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electroesgnetic  pulse.  Paralleling  the  well-known  description  of  laced  circuits] 
In  terns  of  tlvir  poles  and  zeros,  a cocpact  representation  of  the  loop  antenna 
In  tenns  of  its  poles  and  zeros  is  derived.  The  resulting  ti«e  domain  descrie- 
ticn  of  the  loop  response  is  sioply  a suai  of  tene  involving  the  resicurs,  the 
excitation,  and  the  exponentially  Closed  sinusoids  whose  complex  frecuencies  are 
the  frequencies  at  which  red?  alien  fro*  a scatterer  or  antenna  can  take  olace 
without  an  applied  excitation.  One  objective  of  this  research  is  to  investigate 
the  possible  use  of  the  singularity  expansion  method  to  synthesize  radiated  time 
dswiin  waveferws  by  uviforariy  loading  a Iccp  antenna.  In  particular,  tne  wishes 
to  choose  the  loading  si  as  to  realize  some  desired  pole-zero  configuration  on 
the  structure.  Mt  Is  shown  that  the  effect  of  the  loading  can  be  interpreted  as 
introducing  a feedback  loco  into  a block  diagram  representation  of  the  irsedance. 
..transfer  function.  This  observation  permits  one  to  use  the  roct-iccus  techniques] 
} well-known  in  the  area  of  fee  dark  control  theory  to  predict  certain  features  of 
the  pole  trajectories  as  the  leading  ’$  continuously  varied.  Furthermore,  the 
pole  positions  for  a given  impedance  leading  cm  be  found  with  the  aid  of  contdut| 
plots  of  the  Magnitude  and  phase  of  the  impedance  transfer  functie*.  Ccnftlning 
the  use  of  the  above  techniques  for  the  analysis  of  loading  together  with  the 
singularity  expansion  representation,  us  extend  to  electrcraenetie  problems  a 
capability  to  possibly  synthesize  the  desired  response  when  the  input  or  excite- 
jtion  wevefo?*  is  given. 
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A?  though  there  are  many  elertrossacRetic  prebioass 
w&lefc  can  be  treated  satisfactorily  by  weens  of  a static 
or  steady-stats  approximation,  there  are  an  increasing 
auaber  of  problems  in  which  the  transient  behavior  is  of 
p&rsmmt  isportarce.  These  problems  are  usually  diffi- 
cult because  they  pose  the  problem  of  solving  the  field 
equations  as  functions  of  both  tire  and  space. 

A fow  electromagnetic  scattering  and  radiation  prob 
lews  can  be  analytically  solved  directly  in  the  tine 
domain.  However,  for  most  problems  a direct  time  domain 
solution  generally  must  be  obtained  by  numerical  methods 
These  methods  are,  at  best,  tedious  to  apply  and  are 
often  plagued  by  stability  problems. 

A common 1y  more  fruitful  approach  to  obtaining 
transient  field  solutions  is  tc  first  transform  the  tic* 
out  of  the  field  equations.  Host  engineers  art  familiar 
with  this  transform  technique.  In  this  method  the  time 
dependence  Is  transformed  out  of  the  field  equations  by 
either  a Fourier  or  laplare  transform.  The  transferred 
equations  are  functions  of  spare,  with  frequency  appear- 
ing Merely  as  a parameter  of  tne  problem.  The  problem 


Is  then  solved  In  the  frequency  dots* in  either  analyti- 
cally or  numerically  using,  for  example,  a moment  method 
technique.  Once  this  steady-state  solution  of  a problem 
has  been  obtained,  it  is  then  relatively  simple  to  obtain 
the  more  general  solution  representing  the  response  of 
the  object  to  an  Impressed  field  varying  arbitrarily  with 
time.  This  is  done  by  Fourier  inversion  of  the  spectrum 
of  the  solution  quantity  weighted  by  the  spectrum  of  the 
excitation. 

For  solutions  obtained  either  by  time  harmonic  anal- 
ysis coupled  with  Fourier  inversion  or  by  direct  time 
domain  techniques,  a change  in  the  spatial  or  temporal 
behavior  of  the  excitation  requires  that  considerable 
effort  be  spent  in  recalculating  the  response  ©^  the 
structure.  One  is  lead  to  ask  whether  or  not  the  l©* 
established  description  of  lumped  circuits  in  terms  of 
their  poles  and  zeros  might  also  be  used  to  provide  a 
more  compact  representation  of  electromagnetic  field 
problems.  In  the  case  of  electrical  networks,  specifying 
the  finite  number  of  pole  and  zero  frequencies  of  a net* 
work  quantity  (impedance,  transfer  function,  etc.  I com- 
pletely determines  the  quantity  at  ail  frequencies. 
Furthermore,  the  time  domain  response  of  « linear  circuit 
excited  by  an  arbitrary  waveform  may  be  determined  from. 


knowledge  of  the  location  of  these  ride  singularities  of 
the  response  function  in  the  complex  frequency  plane,  as 
well  as  their  corresponding  residues.  The  resulting  tine 
domain  description  of  the  circuit  response  is  simply  a 
sun  of  terns  involving  the  residues,  the  excitation,  and 
the  exponentially  damped  sinusoids  whose  complex  fre- 
quencies are  the  pole  frequencies. 

The  techniques  of  circuit  theory  are  based  on  the 
assumptions  that  path  lengths  in  the  circuit  are  neglible 
and  that  all  electric  and  magnetic  fields  are  essentially 
confined  to  the  circuit  elements.  Field  theory,  on  the 
other  hand,  must  deal  with  fluxes  in  two  or  three  space 
dimensions.  Given  that  circuit  theory  actually  has  Its 
foundations  in  field  theory,  one  night  suspect  that  cir- 
cuit theory  techniques  should  have  analogs  in  the  field 
theory. 

That  the  pole-taro  techniques  of  lumped  circuit 
theory  can  indeed  be  extended  to  electromagnetic  scatter- 
ing was  recognized  by  C.  E.  hau*  [1]  who  for#* Used  the 
singularity  expansion  method  (SEN)  as  applied  to  genera'' 
scattering  problems.  In  his  approach  a conducting 
scatterer  is  described  in  tents  of  an  Integral  equation 
for  «he  induced  surface  current  density.  The  inverse  of 
the  Integral  operator  is  then  expanded  In  terms  of  its 
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poles  end  their  ope r* tor- valued  residues.  The  circuit 
equivalent  of  this  approech  is  the  expansion  of  the 
inverse  of  the  iepedance  matrix  of  an  n~port  network  Into 
a partial  fraction  representation  in  teres  of  the  poles 
of  the  network  and  their  matrix-valued  residues.  Thus  in 

SEM,  field  theory  1$  no  longer  considered  to  bo  something 
apart  free  circuit  and  transmission  line  theories,  but 
rather  as  extensions  of  these  concepts.  Quantities  which 
east  be  known  for  the  expansion  of  the  scattering  operator 
(i  the  inverse  of  the  integral  operator  relating 
induced  currents  to  scattered  fields}  in  terms  of  its 
singularities  art  the  natural  frequencies,  modes,  and 
coupling  coefficients. 

The  natural  frequencies  are  the  frequencies  at  which 
radiation  from  « scattered  or  antenna  can  take  place  with- 
out ar  applied  excitation.  In  other  words,  the  natural 

frequencies  are  the  poles  of  the  structure,  lie  see  lasted* 
otely  that  the  poles  must  be  either  in  the  left  half  o* 
the  s plane  or  on  the  imaginary  axis  in  order  to  exclude 
fields  which  grow  exponentially  with  time,  Poles  on  the 
imaginary  axis,  however,  correspond  to  undamped  sinusoids 
which  therefore  cannot  lose  energy  b*  radiation.  Hence, 
poles  on  the  Imaginary  axis  of  the  $ plane  must  corre- 
spond to  interior  cavity  resonances  which  do  not  radiate 
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exterior  fields.  The  usual  sinusoidal  steady-state 
resonant  frequencies  of  the  structure  are  approximately 
the  leapinary  parts  of  the  coeplex  pole  frequencies. 

For  certain  response  quantifies,  it  is  possible  for  pole- 
zero  cancellation  to  occur.  In  these  cases,  the  natural 
frequencies  do  not  appear  explicitly  in  the  response 
functions  of  the  antenna. 

At  each  pole  frequency  there  is  an  associated  nodal 
current  distribution.  Generally  speaking,  as  a complex 
excitation  frequency  approaches  a natural  resonance  fre- 
quency, the  current  distribution  approaches  that  of  the 
nodal  current  distribution  associated  with  the  pole.  One 
1$  familiar  with  this  behavior  In,  say,  dipole  antennas 
where  at  resonance  the  current  distribution  Is  approxi- 
mately a sinusoidal  standing  wave  with  the  number  of  nodes 
appropriate  to  the  electrical  length  of  the  antenna.  The 
amplitude  of  the  current  depends  on  the  difference  in  the 


pole  and  excitation  frequencies  as  well  as  on  a coupling 
coefficient  which  relates  the  excitation  to  the  proportion 


of  a given  mode  which  is  excited. 

The  objective  of  this  research  is  to  Investigate  the 
possible  use  of  the  singularity  expansion  method  to 
synthesise  radiated  time  domain  waveforms  by  uniformly 
loading  a loop  antenna.  In  particular,  we  wish  to  choose 
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tf»c  lotting  so  « s to  realize  sow  desired  pole-zero  con- 
figuration on  the  structure.  Since  synthesis  design  Is 
usually  carried  out  In  practice  by  Iterated  analyses,  we 
approach  the  synthesis  problem  by  first  building  up  an 
array  of  tools  for  analysis.  These  Include  a rather 
extensive  set  of  tables  for  the  poles  and  residues  of 
unloaded  loops.  The  data  in  thes«  tables  per.it  one  so 

calculate  either  the  tine  domain  or  the  frequency  domain 
response  of  a loop  over  a large  frequency  range  for  an 
arbitrary  excitation.  A product  expansion  representation 
of  the  loop  “transfer  admittance*  function  Is  then 
derived  which  permits  the  rapid  calculation  of  magnitude 
and  phase  contours  for  the  transfer  admittance.  Plots 
of  these  contours,  in  turn,  yield  information  on  the 
shifting  of  pole*  *het  is  possible  by  impedance  loading. 
Adding  further  insight  into  the  problem  of  determining 
the  pole  shifts  are  extensions  of  the  root  locus  tech- 
niques commonly  used  in  control  theory.  The  extensions 
permit  the  techniques  to  be  used  in  the  present  problem 
in  which  there  art  a countably  infinite  number  of  poles. 

With  the  combined  use  of  the  above  techniques,  some 
progress  is  made  toward  the  development  of  an  approach  to 
the  synthesis  problem.  As  in  circuit  theory,  the  Synthes's 
procedure  may  begin  with  either  of  two  different  starting 
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points.  In  the  first,  the  synthesis  problem  is  considered 
solved  when  the  transfer  function  relating  the  response 
quantity  to  the  excitation  has  specified  poles  and  possibly 
specified  residues.  In  the  case  of  the  loaded  loop,  this 
becomes  a problem  of  requiring  the  loading  1 ape da nee 
function  to  interpolate  the  unloaded  loop  transfer  imped* 
ance  function  at  the  pole  frequencies.  If  the  residues 
are  left  unspecified,  it  is  also  possible  to  determine 
whether  or  not  the  synthesised  loading  function  Is  positive 
real , 

The  other  starting  point  sets  out  to  solve  a wore 
difficult  but  wore  practical  problem.  Here  one  is  given 
the  time  domain  response  and  excitation  waveforms  and 
ashed  to  synthesise  the  loading  function  required  to 
approximately  achieve  the  desired  time  domain  response. 

In  this  case*  the  poles  of  the  resulting  structure  way 
not  even  be  needed,  depending  on  the  synthesis  algorithm. 

It  is  eophosirwd  that  the  electromagnetic  synthesis  prob- 
lew  has  an  additional  complication  which  does  not  have  act 
analog  In  lumped  circuit  synthesis.  This  is,  of  corse* 
the  tine  deity  associated  with  the  geometry  c*  the  strut* 
t This  problem  Is  beyond  the  scope  of  this  work  and 
i;  » s anticipated  that  further  development  along  these 
*.1wes  will  rtgui *e  soee  approximation  of  an  infinite 
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CHAPTER  I! 


SINGULARITY  EXPANSION  METHOD  ANALYSIS 
9P  THE  UNLOADED  LOOP 

Tkcrc  exists  8 CSh  tlntilsg  need  to  handle  electro-* 
magnetic  transient  problems  efficiently  and  In  such  a way 
that  a t*fd»  range  of  responses  to  differing  Inputs  may  be 
considered  or  that  the  desired  response  nay  be  synthesized. 
The  purpose  of  this  discussion  Is  to  show  how  the  singu- 
larity expansion  net hod  (SEN)  nay  be  used  to  significantly  . 
simplify  the  calculation  and  synthesis  of  the  response 
of  a transmitting  loop  antenna  excited  by  an  electromag- 
netic pulse.  The  frequency  domain  response  of  a loop  has 
been  extensively  treated  in  the  literature  [2],  [3].  [4], 
aeu  good  summaries  of  these  treatments,  with  some  exten- 
sions* «re  given  by  King  and  Harrison  [S3  and  King  [6]. 

In  the  following,  we  merely  summarise  the  theory  of 
Hu  [4J  as  given  by  King  and  Harrison  [53.  Referral  is 
m tide  to  the  latter  for  details  of  the  derivations,  end 
their  notation  is  generally  followed.  He  have  appropri- 
ately extended  the  theory  of  Hu  [A 3 Into  the  complex  fre- 
quency or,  equivalently,  the  Laplace  transform  domain. 
Although  these  results  may  in,  ohUShed  merely  by  the 
substitution  s » jw  in  Hu's  equations,  we  present  below 
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the  derivation  for  reference  purposes. 

The  usual  method  of  approaching  the  problem  is  to 
write  an  integral  equation  for  the  current  induced  in  the 
loop  which  Involves  the  driving  voltage  waveform.  Because 
sf  the  rotational  *?s2ct ry  of  the  loop,  Fourier  analysis 
of  both  the  excitation  and  the  current  permits  us  to 
derive  a "transfer  impedance"  relating  these  Fourier 
components.  The  modal  transfer  impedance  is  just  the 
ratio  of  the  corresponding  Fourier  components  of  excita- 
tion (voltage)  and  current.  These  transfer  Impedances 
contain  both  the  frequency  and  geometrical  dependences 
of  the  loop. 

2.1  Su»oat*v  of  Mu’s  Theory  for  an  Unloaded  Uecp 

Extended  to  Complex  Frequencies 

I«  the  following,  the  derivation  of  the  solution  for 
the  current  on  a conducting  loop  antenna  I*  summarised 
following  closely  the  presentation  of  King  and  Harrison 
[5],  As  shown  in  Figure  1,  the  center  of  the  loop  coin- 
cides with  the  origin  cf  a cylindrical  coordinate  svster* 
denoted  p,  a,  and  z*  with  the  plane  of  the  loco  lying  in 
the  plane  z * 0.  The  radius  of  the  loop  b is  assumed 
much  larger  than  the  wire  radios  a.  Furthermore,  the 
value  of  a is  seall  compared  with  the  wavelength,  i.e.* 

a*  « b3.  |fca| 2 « 1 (2.1) 
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Tfc«  value  of  !(#')  is  obtained  fro«  an  integral  equation 
obtained  by  invoking  the  boundary  conditions.  The  appro- 
priate boundary  condition  is  that  the  tangential  electric 
field  aust  vanish  at  the  surface  of  the  loop.  The  value 
of  the  electric  field  impressed  across  a delta  gap 
generator  located  at  * - 0 is  ¥*(s).  If  the  structure  is 
i ape dance  loaded,  the  sue  of  the  voltage  drops  across  the 
inpedance  wust  a iso  be  included.  For  a unifornly  loaded 
structure  this  is  easily  acceeplishcd,  since  the  voltage 
drop  per  unit  length  {i.e.t  the  electric  field)  is  werely 
proportional  to  the  current  at  the  same  point.  Consider- 
ing only  the  unloaded  case  here*  we  heve 


v*(sMU) 

b 


(l  f*  * ***) 


» e ■ b 


n.z) 


on  the  surface  of  the  wire  at  #»  where  the  scalar  and 
vector  potentials  at  the  element  d£  « bd$  are  given  by 


$ ■* 


4vc 
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<nV)  »U  * a*) 


(2.3) 
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e')  COS  U - o’)  do'  (2.4) 


Denoting  the  speed  of  light  by  c.  the  kernel  is  defined  by 


Integrating  Equation  {2.1!)  fey  parts* 


il  * — ll-  iL  f l(fe‘)  W{?  - **)  dfe'  (2.12) 
H 4*e$b  H*  J 


and  substituting  (2.12)  and  (2.4)  into  (2.2)*  one  obtains 


«!(S)5(*)  ■ 


J ' 


Mfe  - S')  Ur)  di!  (2.13) 


mfeere  i new  kernel  K{$  - $*)  is  defined  is 

- $')  * -4  COS  ($-$*}-£-.  L-i  lij,*  * 

[ C sb 


U.14) 


kvto  •“$')»  defined  by  (2.5),  e*y  fee  expended  in  a Fourier 

seriws 


W(#  - t')  * }J  k (s)  e 


U.  m 


The  R*  nay  fee  evaluated  in  t*ms  <**  w{#  - $').  yte 


K„(s)  « 


w(*  - *')  ,)j  . * (s) 


Rising  tH«  results  of  {2,15)  and  rewriting  (2.14),  a 
simple  expression  is  obtained  for  K(;  «■  *'!  where 


*-“>  * *=«.**>  * ; |>  (t)  '•  (T) + C«J 

- \ | jV*>  * " 

where  i * «j2sb/e  «nd  I0  and  are  codified  Bessel 
functions  of  the  first  and  second  kind,  respectively,  and 
Y Is  Euler's  constant.  The  constant  Cn  Is  defined  as 


Cn  * In  4n  ♦ y - 2 t 1 

0 Caw  ♦ 1) 


(2,21) 


Using  the  above  results,  (2.13)  reduces  to 


t S<(>  I <(«') 

■*tsr  i 


(Upending  the  above  current  In  a Fourier  series, 


- £ J„U)  e 


iza: 


!n(»)  coefficients  are 


!,(s)  • 


")  O' 


Combining  the  results  of  CS.22)  and  (2.24),  se  obtain 


ijuww 


2 2 *|t*l  (2 .25 


TNI*  1*  * Fourier  saHe*  with  the  coefficients  (joe/2) 
*«C«y  IaU).  The  coefficient*  are  obtained  by  using  the 
properties  of  the  delta  function. 


jne*fiCf)Irt(s) 


it 


L. 
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vj(s)  «U)  t*-*  d* 
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(2.Z6) 


Finally, 


y*) 


n^t^Cs) 


(2.27) 


whore  the  coefficient  -i/Cn^U)]  of  vj(s)  eey  be 
identified  as  the  transfer  ednittance  of  the  nt&  Fourier 
Cttspctteot  of  torrent.  In  the  following  we  repeatedly 
refer  to  this  as  the  "t  master  function"  or 

transfer  admittance.  it*  reciprocal  is  celled  the 


•transfer  iepedance.** 


Z.Z  Expiation  of  the  Transfer  Admittance  in  Terns 

of  Its  Singularities 

To  calculate  the  current,  or  equivalently,  the  trans- 
fer impedance  of  a loop  antenna,  we  have  had  to  solve  the 
corresponding  field  problem-- that  is,  we  have  solved 
Maxwell *s  equations  subject  to  the  boundary  conditions  at 
the  surface  of  the  antenna.  Having  thus  obtained  the 
transfer  admittance  for  the  loop  antenna,  we  next  study 
its  properties  to  determine  the  corresponding  fe^stc 
properties  ef  the  loop  antenna.  However,  some  of  the  west 
important  general  properties  of  the  antenna  transfer 
admittance  may  be  obtained  free  *uo  more  basic  considera- 
tions. These  properties  are  cento i to  all  dynamical 
system-mechanical  and  acoustical  as  well  as  electrical 
sf stems*  and  they  are  independent  of  the  particular  fore 
of  the  equations  as  long  as  these  equations  are  linear. 

Such  properties  were  considered  by  Brune  (7]  with 
special  reference  to  electric  networks,  but  these  results 
are  easily  extended  to  all  linear  dynamical  systems 
Including  systems  with  an  Infinite  number  of  degrees  of 
freedom. 

The  natural  oscillation  constants  of  any  passive 
physical  system,  that  is,  a system  without  concealed 
sources  of  power,  must  lie  either  in  the  left  hatf  of 


the  cOafltx  t ^«m  er  ca  tl«  imagleary  lilts  otl«r«Ht, 
ttt  Ntl  part  of  the  temples  frequencies  of  oscillation 
oeutd  lo  pcs  ft  He,  and  fit  e*ef  Metises  would  grow  fa 
amplitude  without  any  contribution  of  po*er  by  the  system. 
An  antenna  fo  fra*  apace  loans  power  by  radiation, 
whether  ft*  terminals  art  short  efreofttd  or  loft  floating 
banco » the  polos  of  Its  transfer  admittance  am  In  tbo 
loft  half  of  tbo  s piano.  The  only  exception  Is  tbo  point 
at  tbo  origin.  This  point  corresponds  to  a static  field. 

gacegnlzlng  that  the  solutions  of  electromagnetic 

* 

problems  ore  analytic  functions  of  the  couples  free can cy 
s except  at  those  pole  singularities  (couples  natural 
'frequencies  of  oscillation)  Is  the  basis  of  tbo  singularity 
expansion  method  ($£«)  Introduced  by  Sane  £ 1 3 . by 
expanding  tbo  transfer  admittance  In  a partial  friction 
sorlos*  one  needs  only  tbo  polos  and  their  residues  to 
. completely  determine  the  transfer  function  either  in  tbo 
tine  or  tbo  frequency  domain. 

Cue  of  theadvanteges  of  the  singularity  expansion 
method  a*  compared  -to  otbor  more  conventional  methods  Is 
'that  it  provides  « weans  of  characterising  the  electro* 
magnetic  properties  of  a body  with  a discrete  set  of 
Complex'  oumlere  together  with  a set  of  modal  current 
dlsfHbotlobs.  these  quantities  am  uniquely  determined 
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fey  the  bidy  itself  tad  do  act  depend,  for  example,  on  the 
driving  source.  Once  these  quantities  are  fcncwn,  a wide 
variety  of  antenna  problems  can  be  solved  without  having 
to  rt-solve  tuc  boundary  value  problem.  The  singularity 
expansion  method  is  therefore  useful  for  two  reasons: 

(1)  it  provides  physical  insight  into  the  problem  and 

(2)  it  reduces  an  electromagnetics  problem  to  the  minimum 
number  of  quantities  necessary  to  completely  represent 
it. 

The  analytic  property  of  l/an(s)  with  respect  to  s 
allows  the  use  of  various  theorems  of  complex  variables 
in  obtaining  information  about  Its  properties.  The  basic 
Idea  involved  in  this  technique  is  to  expand  the  transfer 
function  of  Van(s)  in  terms  of  its  singularities  in  the 
complex  frequency  plane.  Such  singularities  can  ta&e 
various  forms  such  as  poles,  branch  points  (and  associated 
branch  cuts),  essential  singularities,  and  singularities 
at  Infinity.  For  a restricted  class  of  objects,  which 
includes  the  loop  antenna,  these  s -plane  singularities  are 
limited  to  poles  and  possible  singularities  at  infinity. 

Cnee  one  has  found  the  complex  natural  frequencies 
of  oscillation  and  their  corresponding  current  distribu- 
tions, it  then  remains  only  to  determine  to  what  extent 
each  modal  current  is  excited  by  a given  Fourier  component 
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of  the  input  Ktytf ©r*.  For  the  loop  antenna  the  excita- 
tion of  each  nodal  current  1$  proportional  to  the  product 
of  the  residue  of  the  transfer  adnlttsnce  and  the  Laplace 
transfer*  of  the  corresponding  Fourier  component  of  the 
excitation.  To  deteraine  the  natural  frequencies  of  the 
loop,  it  is  necessary  to  find  the  poles  sn{  of  the  trans- 
fer admittance  factor  1/an(s). 

To  do  this,  one  observes  in  (2.18}  and  (2.20)  that 
an(s)  is  analytic  for  all  s for  n • 0,  and  for  all  s 
except  s * 0 when  n f 0.  This  lapties  that  1/an($)  is 
analytic  for  all  $ except  possibly  for  poles  at  the  zeros 
of  a^U).  Therefore,  since  «n(s)  vanishes  at  these  zeros 
say  Sgf,  Irt(snj}  aay  be  nonzero  with  vn(sn1)  - 8,  i.e., 
no  excitation  1$  required  at  the  natural  frequencies  in 
order  to  have  a current.  Thus,  at  the  pole  frequencies, 
we  have  source  free  solutions  of  the  integral  equation. 
Uaashanhar  fSl  has  shown  that 


*«(s) 


(2.2S) 


$ - s* 


i.e.,  the  transfer  admittance  can  be  written  entirely  as 
a residua  series  involving  its  poles  and  residues  Rni 
This  1$  the  deaf red  expansion  of  the  transfer  ednittance 
in  terns  of  its  singularities. 


By  our  definition  of  the  (tale  frequencies,  ue  hive 
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a*>d  the  corresponding  residues  are  given  by 


Froa  1 2 20)  the  CL«  can  bo  written  as 
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where 
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Using  (2. 26)  in  (2.27),  w my  erlta  the  Fourier 


component  of  the  current  as 
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the  equivalent  expression  in  the  tiwe  domain  Is 


where 


(2.3$) 


and  R^  tad  wist  appear  in  coeptex  conjugate  pairs  in 
order  for  the  tim  doaafn  res  pence  to  be  real.  The  star 
in  (2.34)  denotes  convolution. 

Expression  (2.34)  lap lies  that  poles  oust  be  in  the 
left  half  plane  to  avoid  an  exoonentfally  increasing 
current  as  a function  of  tine.  This  is  explained  phys- 
ically by  the  fact  that  source  free  currents  oust  even- 
tually radiate  away  all  their  energy ; hence 

ift(t)  '£^0,  n M i i0(t)  constant  (2.3b) 
Making  jse  of  (2.23)  and  2.27), 
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(2.37) 


Expanding  aft(s)  in  its  partial  fraction  representation,  we 
have  finally 
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(Z . 38) 


Equation  (2.38)  1$  for  the  c*se  of  an  antenna  excited  at 
4 * 0 by  a delta  gap.  In  the  sere  general  case,  the 
excitation  can  be  represented  by  an  arbitrary  incident 
field  Ejnc($Bs).  Fourier  expansion  of  E*nc(dBs)  aliens 
the  derivation  to  proceed  as  for  the  case  of  an  antenna 
and  the  result  is 
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the  terms  e^  are  called  by  Baum  [1]  the  modal  current 
distributions,  and  the  terms  are  the  coupling 

vectors.  The  quantity 
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is  called  the  coupling  coefficient  and  indicates  hew  much 
of  each  mode  the  incident  field  excites.  The  $„4  are  the 


poles  or  complex  natural  resonant  frequencies. 

If  vj  Is  equal  to  unit  for  frequencies  s * J» 
in  (2. 3?),  then  !(#.j&4  is  the  current  response  due  to  a 
unit  voltage  source  in  the  real  frequency  domain.  Hence 
l(0,jv)  is  the  input  admittance  and  1/If0»j*»)  is  the  input 
Impedance  of  the  antenna.  Hence,  if  a particular  pole 
is  close  to  the  Imaginary  axis*  the  impedance  at  real 
frequencies  In  the  vicinity  of  this  zero  is  snail,  and 
we  have  the  phenomenon  known  as  resonance.  If  several 
poles  are  near  the  Imaginary  axis,  the  lupedance  will 
fluctuate  between  i«all  and  large  values  as  the  frequency 
passes  these  points  As  the  poles  recede  from  the  imagi- 
nary axis,  the  fluctuations  beeves*  less  pronounced,  and 
the  “resonance  curves"  become  flatter. 

An. equivalent  circuit  can  be  developed  representing 
the  field  equations  of  Kaxwell  for  an  electromagnetic 
field  containing  conductors  and  bound  charges.  Both 
transient  and  sinusoidal  field  phenomena  may  thus  be 
studied  by  numerical  and  analytical  circuit  methods; 
such  « circuit  modal  applies,  of  course,  to  radiation 
from  a loop  antenna  [9].  The  circuit  models  can  in 
principle  bo  developed  for  all  curvilinear-orthogonal 
reference  frames  to  allow  the  solution,  to  any  desired 
degree  of  accuracy,  of  any  two  or  three-dimensional 


problem.  The  models  correspond  to  the  approximation  of  e 
transmission  line  In  one  dimension  by  • casceded  series  of 
sections  containing  ordinary  luaped  circuit  eleaents  R,  L, 
C,  and  G's.  Since  the  field  equations  of  Kaxwell  aay  thus 
be  represented  by  a stationary  network  (within  any  desired 
degree  of  accuracy),  it  way  be  stated  that  [9]: 

Any  theorem,  formula,  concept,  or  law  that 
is  valid  for  stationary  networks  (such  as  reci- 
procity theorems,  Thevenin's  theorem,  concepts 
of  dualism,  reduction  formulas,  generalization 
postulates,  etc.)  can  be  translated  into  a core- 
sponding  theorem,  formula,  concept,  or  law 
relating  to  the  electromagnetic  field. 

As  the  number  of  elements  in  the  circuit  increases, 
the  number  of  poles  also  increases.  As  3 consequence  of 
the  above  circuit  model  of  Maxwell's  equations,  contin- 
uous structures,  including  all  of  free  space,  are  * twits 
of  networks  with  an  increasingly  larger  number  of 
increasingly  smaller  meshes.  The  number  of  their  zeros 
and  poles  will  be  infinite.  {In  fact,  this  must  bo  true 
of  any  physical  circuit  since  all  physical  circuits  **■« 
continuous  and  cannot  be  wholly  disassociated  from  the 
surrounding  space.) 

Thus,  there  are  an  infinite  number  of  complex 
resonant  frequencies  whose  locations  are  in  the  left  half 
of  the  complex  frequency  piano,  and  which  occur  in  complex 
conjugate  pair'  as  a consequence  of  the  above  theorem. 


la  order  to  generate  a real -time -domain  response,  pdas  of 
the  transfer  impedance  function  must  appear  in  complex 
conjugate  pairs  and  their  residues  must  be  in  complex 
conjugate  pairs  for  poles  not  on  the  negative  real  axis. 

It  is  then  recognized  that  the  transient  response 
of  an  object  can  be  viewed  as  a superposition  of  a series 
of  deeped  sinusoidal  oscillations  at  the  so-called  natural 
frequencies  of  the  object,  it  has  been  observed  in  cany 
electromagnetic  pulse  (EKP)  scattering  and  interaction 
problems  that  the  time  dependence  of  various  quantities* 
Such  as  the  current  induced  on  an  object  in  an  ENP  simu- 
lator, seems  to  be  described  by  only  a few  cf  these 
exponentially  damped  sinusoidal  csci liationsCl]. 

2.3  huserical  Techniques  and  Results 

The  expression  for  e^s),  (2.16)  to  <2,21),  involves 
integrals  of  Anger-lieber  functions  of  complex  arguments. 
These  functions  have  been  computed  using  an  extension 
to  complex  arguments  of  the  methods  described  in  £10]. 

A parametric  study  of  the  roots  of  att{x}  as  a function 
Of  the  ratio  b/a  has  been  carried  out  using  a numerical 
search  procedure.  The  roots,  of  course,  satisfy  the 
requirement  that  they  should  appear  only  in  the  left  half 
of  the  s plane  and  in  complex  conjugate  pairs. 

The  so-called  (fuller's  method  has  been  used  to 


Mjaerieal ly  determine  the  roots  of  an(s}.  In  this  method 
whret  values  cf  an(s)  are  computed  from  three  estimates 
of  the  root,  say  at  sn-2,  sn_i»  and  sR.  Using  these 
values,  a quadratic  interpolation  formula  [H]  is  used 
to  approximate  »n(s)  in  the  vicinity  of  the  given  points. 
The  root  of  the  quadratic  nearest  the  best  estimate,  say 
*n.  is  designated  and  the  procedure  is  repeated  using 

Vl*  V sn*l  until  ^$a*l  * V / ^n  * U ]<6%$  th*n  a 

preaasigned  nuRber.  For  more  details,  the  reader  is 

referred  to  reference  [11]. 

Approximate  pole  locations  to  initialise  the  proced- 
ure are  net  necessary,  though  in  practice  having  good 
initial  estieates  of  pole  locations  will  considerably 
reduce  the  number  of  iterations  required  to  find  the 
various  roots.  This  searching  procedure  Is  in  general 
found  to  5c  quite  efficient,  but  it  has  bees  found  to  be 

rather  difficult  to  find  all  of  the  poles  in  certain 

regions  when  good  initial  estimates  ara  unavailable. 

Ore  method  which  has  been  employed  to  overcome  th?s 
difficulty  H to  actually  plot  magnitude  contours  of  the 
function  in  the  complex  s plane.  The  region*  near  ter  os 
of  an(s)  Shew  up  clearly  on  contour  plot  and  provide 

more  accurate  initialising  cat*  for  Muller’*  Method. 

Furthermore,  the  contour  plots  themselves  are  useful  in 
determining  the  poles  of  the  loaded  structure,  as  *<*  sha^ 


wv 


mtmiu 


i@©  in  Chapter  111. 

The  residue  for  1/a^Cs)  m$  be  easily  calculated  by 
the  residue  theore*  using  a circular  coo tour  about  the 
pole  as  shove  In  Figure  Z . 

By  the  residue  theorem,  we  have 


fi 


ds 


2»j  *n(s) 


C2.41 


la  Figure  2,  let 


s - S|  * ee^*. 


s * 5^  ♦ ee* 


ds  * jeeJ^  do 


(2*42 


and  substituting  (2.42)  Into  (2.41),  we  have 


jce^  do 


a a — — up  -«<r .■■—-*.»»- 

2*J  Jcnian(s)(sn1  ^ ce 


For  a nu«@rica1  approximation,  we  say  divide  the  contour, 
Cn|  into  « equal  subdivisions  and  apply  the  s tuple  reef, an 
gular  rule  for  integration: 


m 


“ f-;  E “ 


j 2(n-l)t/» 


2*  n~  *rt(s)  ♦ ce*  2(n-1)e/s 


(2.44) 


Thus  the  residue  at  the  pole  zni  for  a (s)  is 


(2.45) 


e ■ £(n-1)v/n 

* ■ S «„<»)[.,  ♦ *C»-1  )«/-J 

2s  equation  (2.28)  an  important  feature  Is  that  the  func- 
tion l/aaCs)  may  be  represented  by  Us  poles  and  residues 
with  no  additional  entire  function  required,  as  shown  by 
Uft&ihaakar  [8]. 

Several  checks  were  mate  to  determine  the  rate  of 
convergence  of  (2.4S).  The  value  of  e «as  varied  fro* 
!<Tl  to  10“s  while  simultaneously  the  value  of  a was 
varied  frm  3 to  24.  it  was  found  that  for  the  combina- 
tion of  c * 10*  * and  ■ • 3,  accurate  answers  with  a 
relative  error  of  the  order  of  10*’  were  obtained.  4 
large  number  ef  these  pole  locations  s * $e|  and  corre- 
sponding residues  8^  Have  been  calculated  and  tabulated 
for  a wide  range  of  the  loop  parameter  8 * ZZn(Ub/a}. 
These  results  permit  one  to  use  (2.34)  to  calculate  time 
domain  loop  currents  for  arbitrary  excitation  without 
resorting  to  the  comparatively  inefficient  process  of 
Fourier  transformation. 

By  direct  calculation  we  find  in  the  eo&plex  fre- 
quency plena  an  infinite  number  of  camples  resonant 
frequencies  whose  locations  suggest  three  separate  cate- 
gories of  resonant  frequency  or  pole  types  for  each  node 
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yy,-f 


■c 


*%2-% 
6/c,  ■>  • ■ ,| 

K •■ 


■ rr 


Sfcrjl 


i^’  ' ;•;•->>■  ' 


„„J| 

SPfR.i 

*Sv5? 


J$@-- 


'Wfr'-K-i 


v ^ 

r>^PJ5-^**V: 


HSL 

Tfccrt  Is  a single  pale  very  near  the  s * jw  axis  at 
approxisately  w » n,  This  pole  gives  the  principal  con- 
tribution to  the  ties  domain  response  of  the  loop  at  late 
tfsaes  and  the  imaginary  part  of  the  pole  location  corre- 
sponds closely  to  the  resonant  frequency  of  the  loop  for 
an  excitation  of  the  fere  e^8*. 

Type  II 

There  arc  n*1  of  these  poles  ( including  cen jugate 
pairs!  which  lie  roughly  on  the  left-hand  side  of  an 
ellipse  centered  at  i » 0 and  with  e seafssajor  axis  some- 
what larger  than  a. 

Typ*  m 


There  is  a layer  of  poles  lying  alamt  parallel  to 
» s * jy  axis.  The  layer  contains  an  infinite  nueoer 
of  pgles  and  they  a fig  spaced  approxiaately  &m  « *c/b 
units  apart,  where  b is  the  loop  radius, 

JU  with  thin  cylindrical  wires,  increasing  the  wire 
radius  hat  the  effect  of  shifting  the  Type  i poles  near 
the  ju  axis  away  free  the  axis,  or  egui valent 1y«  increas- 
ing the  deeping  constants  of  those  nodes  in  the  tine 
dewain  [12].  Types  £1  and  III  poles  located  further 
eway  front  the  imaginary  axis,  however,  wove  away  from  the 


'M 


t 


imaginary  axis  as  the  radius  approaches  zero. 

Pole  and  residue  data  for  the  three  loop  sizes 
a • 2£i2*b/a  * 1C,  15,  and  20  are  presented  in  Figures  3 
through  8 and  Tables  1 through  6 for  modes  n * 0 through 
40.  For  a particular  mode,  Type  I!  poles  fall  on  an 
elliptically  shaped  curve  with  n *■  1 poles  (deluding 
conjugate  pairs).  There  will  be  one  more  pole  (Type  1* 
at  approximately  a • n.  Displayed  in  Figures  4,  6,  and  8, 
corresponding  to  fi  *•  10,  IS,  and  20,  respectively,  are 
the  layers  of  Type  HI  poles  parallel  to  the  s » jw  axis. 
These  poles  are  shown  for  each  mode  0 through  20  for 
values  of  ub/c  « 0 to  ob/c  * 36. 

The  residue  corresponding  to  each  of  the  poles 
plotted  if*  Figures  3 through  8 are  tabulated  in  Tabl#$  1 
through  8.  These  are  tabulated  in  three  colusaas  far  the 
three  loop  sizes  Q « 16,  IS,  and  20.  The  first  number  in 
each  column  represents  the  real  value  of  the  residue,  and 
the  second  is  its  imaginary  value.  The  Index  is  a imigut 
four-digit  number  that  identifies  & particular  pole.  The 
first  two  digits  are  the  mode  number  for  erodes  0 through 
20.  The  last  two  digits  tan  be  grouped  into  three  differ- 
ent categories  corresponding  to  the  three  different  types 
of  poles.  A double  zero  (00)  in  the  last  two  digits  of 
the  index  corresponds  to  the  single  Type  I pole  very  near 
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Tititf  3.  test  dues  of  the  Poles  fo*“  Pole  Indices  1206  to  1601 
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tluk  akfc  axis  at  approximately  « » a.  Type  II  pelts, 
which  vary  tn  number  according  te  the  soda  number  n,  have 
indices,  tlia  last  tuo  digits  of  which  range  from  01  to  as 
high  as  Tin 

For  example,  for  0 * 10  node  20  (Figure  3),  the 
index  number  2009  represents  the  pole  at  (-1.08,  19.96) ; 
2901  raprasaats  the  pola  at  VI. 77  , 21.82);  2002  the  pole 
at  (-5.70,  16.99);  aad  finally,  2011  the  pole  at  (-13.66, 
0.9).  For  a given  node,  aa  Increasing  nuaber  In  the  last 
ta(€»  index  digits  novas  along  the  elliptical  arc  fro*  near 
the  wb/c  axis  toward  the  negative  real  axis  In  a ecus ter- 
clockwise  manner.  TSies^  first  two  types  of  poles  Have 
Own  plotted  in  Figures  3,  S,  and  1.  The  third  category 

it?  poles.  Type  IS  I,  contains  m Infinite  nuaber  of  poles 
’yinp  almost  parallel  to  the  ub/ c axis.  Only  the  poles 
such  that  wb/c  « 39  are  tabulated.  Again,  the  first  two 
dif  'ts  of  the  1**dex  *w-bo«*  represent  the  mode  number; 
however,  bc.-a  Va  ’ o ; two  *!1gits  all  begin  with  3, 
denoting  p«*  :•$.  of  r-  - ^ type.  These  range  from  31 

up  to  as  *n>%  as  **.  For  example,  for  D * 10  mode  0 
(Figure  4)  toe  i«0£*  leader  0931  represents  a pole  at 
(-0.69,  3,76);  993?  represents  a pole  at  (-9.93,  7.35); 
and  9039  a pole  at  c-l.lS,  29.66).  Poles  in  this  layer 
are  numbered  sequentially  beginning  with  the  pole  having 
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ta«  seal left  imaginary  part  and  proceeding  away  frow  the 
o axis.  These  Type  III  poles  are  plotted  in  Figures  4,  6* 
and  8. 

Uaa Shankar  [8]  notes  that  asymptotically  for  Type  III 
poles  there  are  only  two  sets  of  roots*  one  set  for  even 
codes  and  another  set  for  odd  sodas.  The  convergence  to 
these  two  sets  of  values  for  the  lowe*  orders  can  be 
readily  seen  In  Figures  4.  6,  and  8.  “her*  for  Increasing 
values  of  s along  the  mb/ c axis  all  Type  HI  poles  con- 
verge. regardless  of  mode,  to  one  distinct  set  of  values 
for  even  nodes  and  one  for  odd  nodes.  In  Tables  1 through 
8 it  can  oe  seen  that  for  relatively  large  values  of  $ 
the  residues  follow  the  sane  pattern.  For  large  s all 
even  anodes  tend  to  the  same  residues  regardless  of  node; 
similarly  for  all  odd  codes. 

These  tables  of  pairs  of  cowplex  numbers  representing 
the  polos  and  residues  of  the  admittance  transfer  function 
are  uniquely  determined  by  the  loop  geometry.  Independent 
of  excitation.  They  provide,  through  the  partial  fraction 
expansion  of  the  admittance  transfer  functions,  a means  of 
accurately  characterising  the  electromagnetic  properties 
of  loop  antennas  for  throe  relatively  different  slstd 
antennas  through  mooe  20.  The  representation  of  the  time 
domain  response  in  terms  of  the  poles  t_.  and  residues  SL< 
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yields  time  domain  results  which  coupe  re  veil  with  the 
Fourier  transferal  of  frequency  dos&in  data.  However,  a 
comparison  in  the  frequency  domain  of  the  integral  repre- 
sentation and  the  partial  fraction  representation  of  the 
transfer  admittance  function  at  the  same  frequency  indi- 
cates iitit  the  latter  representation  does  not  appear  to 
converge  to  the  correct  result.  This  is  Illustrated  in 
Figure  9 where  the  partial  fraction  and  integral  repre- 
sentation of  l/a0(js4  are  compared.  One  notes  the 
apparently  constant  offset  in  the  imaginary  part,  which 
should  he  zero  it  u * 0.  Since  a constant  would  repre- 
sent an  additional  entire  function  added  to  the  partial 
fraction  expansion,  this  problem  was  studied  in  some 
detol  1 . 

An  obvious  possible  source  of  error  would  be  that 
the  sn|  and  f?n^  were  not  being  computed  to  syf*1  cle.it 
accuracy,  the  early  results  were  computed  using  standard 
precision  arithmetic  on  the  CDC  6600  with  H significant 
digits  of  accuracy.  The  Anger- Weber  function  calculation 
was  evaluated  by  using  m alternating  series  expansion. 

To  check  ?er  roundoff  errors,  the  entire  routine  was 
rewritten  for  double  precision  which  carries  28  signifi- 
cant digits.  The  new  routine  was  checked  against  the 
original,  but  the  troublesome  imaginary  value  at  the  origin 
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remained  virtually  unchanged. 

One  insult  of  the  conversion  to  double  precision  was 
the  extension  of  the  range  of  the  argument  for  which  the 
Anger-ieber  function  could  be  accurately  calculated.  This 
permitted  poles  sRl  and  residues  Rai  to  be  calculated  In 
the  region  M/c|  < 30.0  in  the  s plane.  Using  the  result 
of  [83,  the  asymptotic  formula  was  also  used  to  compute 
sni  *n<!  for  values  near  db/c  * 30.0.  In  the  regions 
where  the  two  calculations  could  be  compared,  soma  small 
differences  were  noted  between  the  accurate  value  of  sfl| 
and  calculated  using  th*  series  expansion  method  and 
that  computed  using  the  asymptotic  formula.  It  is  cer- 
tainly true,  then,  that  for  values  of  sn^  and  8^  4n  the 
range  of  db/c  from  30  up  to  possibly  several  Hundred, 
small  errors  are  introduced  by  the  asymptotic  formula. 

At  s * 0,  the  partial  fraction  series  becomes 

!/aa(o)  • - E Ra|/sn<  (2.48) 

where  the  series  must  sum  to  zero  to  have  the  correct 
value.  It  is  easily  seen  that  errors  in  the  poles  and 
residues  will  give  l/an(s)  the  wrong  behavior  near  s «*  2. 
Furthermore,  dividing  the  partial  fraction  series  Into 


those  terms  whose  poles  and  residues  are  computed  from 


tat  series  representation  of  the  tager-Veber  function  end 


those  computed  froa  the  esyaptottc  formula,  we  heve 
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The  lest  series  for  the  range  of  s ■ j**,  considered  In 


Figure  9,  can  be  approximately  replaced  by 


f 

’s.*b> 


(3.48) 
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since  s is  small  coopt  rod  to  s„4 . Errors  in  these  terns 


would  explain  the  constant  offset  sheen  in  Figure  9.  Add- 


ing further  weight  to  this  argument  is  the  fact  that 


sittilar  plots  of  l/afi(s}  for  other  values  of  « in  the 
seise  range  of  % shoe  the  sans  nearly  constant  offset » the 


value  of  the  constant  being  independent  of  n.  This  Is  not 


surprising,  however,  when  one  recalls  tft««  tU  c;>*vtctic 
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formula  predicts  the  sice  set  of  poles  and  residues  for 
ill  even  n and  another  set  for  all  odd  n. 

In  order  to  alleviate  the  conversance  difficulty  In 
the  s domain*  an  Infinite  product  representation  was  used 
The  product  representation  has  the  advantage  that  both 
the  poles  and  zeros  of  the  function  are  automatically 
Included  In  the  representation.  A further  advantage  of 
this  formulation  Is  that  the  residues  *nl  ore  not  require 
These  expansions  for  t/an(s}  are  derived  in  Appendix  A 
and  are  given  by  (2.49)  and  (2. SO): 
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for  « t 0 where  the  are  singularities  of  the  n 
cod  is  evaluated  at  s equals  zero; 
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when  « « 0 or  where  are  the  singularities  of 
and  hs  U evaluated  at  s equals  zero. 


mode  zero 
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Another  calculation  of  l/ae(s)  uutag  (2. SO)  appears 
in  Figure  9.  The  offset  In  the  imaginary  part  Is  coa* 
plenty  eliiifnctad.  By  reforeualtlag  the  problem,  the 
constant  offset  has  disappeared,  and  the  agreeaent  of  the 
series  and  product  representations  t£2d  to  validate  the 
accuracy  of  the  calculated  sn{.  It  Is  therefore  concluded 
that  the  poor  convergence  of  (2.28)  In  the  frequency 
do as in  is  due  to  seall  deviations  in  the  computed  location 
of  the  s_4  and  values  of  *L,  for  wh/c  > 39.0  in  the  s 

B1  fn 

plane. 

The  product  expansion  provides  a rapid  and  accurate 
means  of  calculating  the  values  of  l/*n(s)  in  the  complex 
frequency  plane,  whereas  the  Laplace  inverse  of  the 
partial  fraction  expansion  is  useful  and  accurate  in  the 
tine  dosain.  Those  two  representations  provide  the 
necessary  tools  to  aecoaplish  all  required  calculations 
for  the  loop  antenna  accurately  and  quickly. 

Sue  further  question  remains  concerning  the  poles 
for  the  unleaded  loop,  it  has  been  pointed  out  that 
the  third  category  of  poles,  l.e..  the  ones  per* tie  ling 
the  imaginary  axis  in  the-  s plane,  do  not  correspond  ce 

similar  poles  of  either  dipoles  or  spheres.  Seme  specu- 
lation has  been  rede  that  these  poles  are  due  to  the  thin 
wire  approximations  used  to  derive  the  transfer  functions 
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*B(s).  0 «e  should  keep  in  aind,  however,  that  such  struc- 

tures as  dipoles*  spheres,  and  prolate  spheroids  are 
topologically  identical  in  that  one  structure  can  be 
continuously  deforced  into  another  and  that  as  such,  there 
oust  be  a one-to-one  correspondence  of  their  poles.  This 
corre^ondence  is  deterained  by  noting  which  poles  aerge 
Into  those  of  a sphere,  say,  as  an  object  is  continuously 
deforced  into  a sphere.  The  loop,  however  is  net  topo- 
logically equivalent  to  a sphere,  however,  because  a 
sphe/e  without  handles  can  never  be  deforced  into  a Teoo 
and  its  poles  do  not  necessarily  correspond  to  these  of 
a sphere.  To  test  this  hypothesis,  a ssethod  of  moments 
solution  for  the  f*  * 0 mode  current  induced  in  a toroidal 
antenna  of  2 * 10  was  implemented.  The  aathematical 
derivation  of  the  integral  equation  together  with  the 
numerical  considerations  are  contained  in  Appendix  S. 
Basically,  the  loop  is  treated  as  a conductive  toroid 
divided  up  into  a large  rtus&er  of  curvilinear  patches. 

The  wire  circumference  was  d .ided  into  2«  segments  with 
the  current  assumed  to  be  urifer*  in  each  secant.  Since 
in  Mu's  solution  the  function  1/eM(s)  is  proportional  to 
the  total  current  for  a uniformly  excited  wire,'  in  the 
moment  solution  the  current  density  on  the  wire  was 
integrated  to  find  the  total  current  under  the  condition 
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cf  sfsifore  excitation  of  the  satire  wire  surface.  The 
total  current  found  this  way  corresponds  to  l/ac($)  ?nd 
tha  contour  plots  of  t^e  current  should  yield  the  true 
poles  of  the  loop  {within  Benefit  method  approximations}. 

It  was  found  that  in  order  to  generate  the  required 
values  for  filling  the  ifcatrix,  iso  re  computer  tine  was 
required  than  expected.  Consequently*  it  was  necessary 
to  ensoloy  a minimus  number  of  points  in  the  contour 
generation  scheme  in  order  to  keep  the  required  computer 
'Mfc*  within  reasonable  limits.  The  resulting  transfer 
ad«.1  for  node  0,  n * 10  is  shown  in  Figure  10.  ’^e 

will  sSsik  . n Chapter  U!  that  by  comparing  phase  and 
scagnitude  values  near  the  isav  nary  axis  for  the  method 
of  moments  and  the  product  expansion  {Figure  11)*  one 
finds  nearly  the  same  pole  structure.  The  poles  sQi  of 
Type  HI  are  not  only  present,  but  are  in  approxl  ?sately 
the  casre  location  as  predicted.  However,  a set  of  zeros 
in  the  total  current  lies  interspersed  with  the  poles  in 
a zigza-j  fashion  running  paralles  to  the  imaginary  axis. 
This  ',ct  of  zeros  in  the  total  current  Is  caused  by 
Surface  currents  which  flow  in  opposite  directions  on  the 
inside  and  outside  of  the  wire  such  that  the  total  current 
is  zero.  In  addition  to  the  zeros,  there  is  also  a second 
layer  of  pcle-zero  pairs  which  parallels  the  firs*,  ilhi  le 
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this  layer  has  such  large  damping  constants  that  It  has 
negligible  contribution  in  the  tiae  and  frequency  domain, 
its  appearance  is  interesting,  nevertheless. 


The  matrix  deterclnant,  of  course,  would  not  have 
the  zeros,  but  in  order  to  plot  the  determinant,  the 
density  of  points  calculated  would  have  had  to  be  increased 
by  10  to  100  tines.  This  is  because  the  appearance  of 
the  zeros  of  the  determinant  were  found  to  be  extremely 
localized  in  the  s plane. 

ho  attempt  was  made  to  generalize  the  moment  method 
program  to  Investigate  the  Type  III  poles  for  modes 
other  than  n * 0,  but  it  is  expected  that  they,  too,  do 
indeed  exist. 
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CHAPTER  III 


ANALYSIS  OF  THE  UNIFORMLY  LOADED  LOOP  AJflTEH.IA 

The  analysis  of  an  unloaded  loop  an*~nna  Has  ueen 
presented  in  the  previous  chapter  using  the  singularity 
expansion  Method.  in  the  present  chapter,  the  effects  of 
uniferaly  loading  a conducting  loop  are  considered.  Frew 
the  point  of  view  of  the  singularity  expansion  aethod, 
the  loading  aerely  shifts  the  location  of  the  odes  of 
the  structure.  It  is  shown  that  contour  plots  of  the 
transfer  iepedance  defined  in  Chapter  II  Bay  be  used  to 
find  the  snifted  pole  positions  end  to  determine  what 
shifts  are  possible.  Furtherstore,  it  is  shown  that  the 
effect  of  the  loading  can  be  interpreted  as  introducing 
a feedback  loop  into  a tncck  diagraa  representation  of 
the  impedance  transfer  function.  This  observation  permits 
one  to  use  the  root-locus  techniques  well-known  in  the 
are*  of  feedback  controls  tc  predict  certain  features  of 
the  pole  trajectories  as  the  loading  is  continuously  varied. 
Finally,  representative  tire  domain  calculations  for  the 
step  response  of  a linear  antenna  are  given  for  various 
values  of  purely  resistive  loading. 
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3.1  Derivation  of  the  Admittance  Transfer  Function  for 
the  Loaded  Loop 

Consider  uniformly  loading  the  loco  with  an  impedance 
ZL(s)  total  Impedance  around  the  loop  or  Z^(s)/2*b 
impedance  per  unit  length.  The  boundary  condition  on  the 
loop  is  that  the  tangential  electric  field  equal  the 
product  of  the  total  current  and  the  impedance  per  unit 
length.  Thus  the  integral  equation  (2.2)  is  replaced  by 
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Expanding  J,  A , and  l(o)  as  in  section  2.1,  we  have 
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which  can  be  written  as 
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a for*  which  parallels  (2,25).  Hence  corresponding  to 
(2.27),  we  now  have 
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n»  aR(s)  - jZL(s)/ntfen 

This  transfer  function  has  poles  at  frequencies  s e 
such  that 

Msil)  * !3'6) 

or  equivalently, 

j*n(sni)  ! ■ kMH  (3-J) 

■ ‘ 90”  <3-8> 

where  the  prime  distinguishes  quantities  defined  for  the 
loaded  loop  free  those  defined  for  the  unloaded  loop. 

The  argument  given  by  Uwashankar  fS]  for  the  expansion 
of  the  admittance  transfer  function  applies  here  also 
and  results  in  th&  partial  fraction  expansion 
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The  residues  8^  In  (3.9)  above  are  easily  found  to  be 


(3.10) 


where 
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(3.11 
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The  fore  of  ttee  partial  fraction  expansion  1 1 (3.9)  above 
Is  Identical  to  (2. 28}  but  cith  priwed  quantities  replacing 
unpriced  quantities.  Hence  the  discussion  cf  the  tlsae  and 
frequency  domain  current  response  ’n  Chapter  U applies  to 
th*  loaded  loop  at  well. 

3.2  Use  nf  Contour  Plots  to  Represent  Poles  cf  the 
leaded  loop 

Equations  (3.7)  and  (3.8}  indicate  that  contour  plots 
c*  the  magnitude  and  phase  of  a (s)  in  tee  cc*oiex  frequency 
olar.e  would  simultaneously  be  plots  of  the  magnitude  anc 
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phase  (shifted  by  96*)  of  the  normalized  Ispedartce  loading 
required  ef  each  point  In  the  complex  fluency  icnain  In 
order  to  have  a pole  there.  Using  the  product  expansion 
foreulas  (2.49)  and  (2.50),  and  including  all  poles  In  the 
range 

|la^$njb/cj  j <,  200  (3.13) 


Magnitude  and  phase  contours  of  the  afl(s)  are  plotted  in 
figures  11  through  21.  these  figures  are  for  Modes  n * 0 
through  n • 20  for  a * 10.0. 

the  Magnitude  lines  m * constant  in  the  figures  deter- 
*ine  the  contours  of  constant  Magnitude  of  §n  and  the 
nemsallaed  Inge  dance  loading  according  to 
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(3. IS) 


The  phase  o in  degrees  is  given  for  Zt  {%)  and  is  related 

fc* 

to  the  phase  of  a„(s)  by 

0 • argjlt(s)j  • argj@B(s)]  - SO*  (3. IS) 


k n usher  of  observations  concerning  the  contour  plots 
are  appropriate.  One  notes,  for  instance,  that  only  ae(s) 
has  a pole  at  s * 6;  for  all  other  *„(*)•  s • 0 is  the  only 
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zero  of  the  transfer  impedance.  Furthermore,  one  notes 
that  at  each  point  on  the  contour  plot,  we  simply  read  off 
the  magnitude  and  phase  the  impedance  should  have  In  order 
to  have  a pole  at  that  point.  This  fact  Is  very  Important 
and  Is  exploited  In  the  next  chapter  In  synthesizing  the 
loop  response  by  choosing  the  pole  locations.  When  the 
loading  function  2j_ f s ) is  given  (i.e.,  in  the  analysis 
problem),  some  additional  effort  Is  required  to  graphically 
find  the  pole  locations.  Mhat  is  needed  1$  a separate 
contour  plot  of  ZL(s)  using  the  sane  complex  frequency, 
magnitude  and  phase  scales.  Then  by  overlaying  the  two 
plots,  loci  of  common  magnitudes  and  phases  can  be  drawn. 
Their  intersection  will  be  the  pole  positions  for  the  given 
loading  function  Z^fs).  ^ote  that  this  would  have  to  be 
done  for  every  ssode,  » D,  1,2,  . . Decause  the  lead- 
ing affects  all  poles. 

Ore  notes  that  it  is  oossiole  to  shift  poles  into  the 
r - jisi  half  plane,  but  this  could  only  be  done  with  active 
loading.  However,  generally  speaking.  Figures  11  through 
21  indicate  that  the  heavier  the  loading  (that  ss,  the 
larger  the  magnitude  of  lAs))  the  larger  the  damping 
constant  becomes  for  the  shifted  pole.  Cn»  should  recall, 
however,  that  in  Figure  13  In  Chapter  11,  which  was  com- 
puted usiag  jcossmnt  methods  with  ss  thin  w'-e  assumptions. 
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one  finds  in  actuality  & layer  of  zeros  just  beyond  the 
first  layer  of  poles  for  node  n * 0 and  that  heavy  load- 
ing shifts  the  poles  toward  these  zeros  rather  than  toward 
infinity.  A similar  situation  will  probably  exist  for 
higher  order  aodes. 

finally,  it  is  noted  that  the  lines  where  the  phase 
is  zero  on  each  plot  correspond  to  purely  resistive  load- 
ing and  that  if  the  resistance  is  frequency  Independent, 
then  the  poles  will  all  be  on  the  same  magnitude  contour, 
further,  note  that  the  Type  I pole  (at  s ■ 0)  and  the 
Type  II  pole  (on  the  negative  real  axis)  for  mode  n « 0 
approach  each  other  with  increased  resistive  loading, 
fore  a double  pole  on  the  negative  real  axis  for  a certain 
critical  resistance,  and  then  sol  it  into  a complex  con- 
jugate pair  with  increased  loading.  The  dominant  Type  I 
pole  for  eode  n - 1 approaches  the  negative  real  ax: 
with  increased  resistive  loading,  forms  a double  poie 
there  with  its  conjugate  pair,  and  then  the  two  poles 
solit  apart  and  move  away  from  each  other  alone  the  nega- 
tive real  axis.  These  are  tne  only  two  cases  'or  which 
dominant  (Type  I)  poles  combine  with  other  poles  to  form 
a double  pole.  At  either  of  these  loading  conditions, 
one  may  call  the  antenna  *o-1ttcaliy  damped. * The  trajec- 
tories of  the  dominant  poles  for  modes  n ® 3 to  « * 7 
with  resistive  loading  only  are  shown  in  f i dure  22. 


3.3  Feedback  Interpretation  of  Impedance  loading  and 
Root  Locus  Methods 

The  relation  between  the  excitation  and  the  Fourier 

components  of  current  f3«-  the  unloaded  loop* 


-i  Ve(s* 
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aay  be  represented  in  block  diaorae  fore  as  in  Figure  23 
The  corresponding  relation  for  the  loaded  loop  is 


JZ,  {sj/ttm 


which  can  be  represented  In  the  block  diagram  of  Figure  2*. 
One  sees  tnat  the  effect  of  the  loading  is  to  add  a feed- 
back loop  into  the  unloaded  loop  transfer  function  which 
shifts  the  poles  of  the  original  unloaded  system.  This 
interpretation  of  leading  as  adding  «?  feedback  path  oerwits 
one  to  use  the  techniques  of  control  syste*  theory  {12} 
to  analyse  * for  example,  the  effect  of  feedback  on  the  pole 


locations  and  their  sovenent  as  the  feedback  (loading) 
changes,  tie  discuss  in  the  following  the  use  of  the  root 
l«c,*s  technique  of  control  theory  [14]  as  applied  to  the 
los&H  loop. 


In  order  to  apply  the  root  locus  techniques » it  is 
convenient  to  think  of  attaining  a certain  given  uni  fore 
loading  ZL(s)  by  loading  the  loop  with  a uni  fore  ispedanee 
where  the  real  constant  & is  varied  fro*  £ * 0 
(corresponding  to  an  unloaded  loop)  to  £ • 1 (corresponding 
to  loaded  with  the  desired  loading).  To  consider  infinite 
loading,  we  *ereiy  let  £ tend  to  infinity.  Thus  the  loop 
transfer  function  for  a given  loading  £2L(s)  is 


(3.19! 


The  poles  of  the  systew  are  now  those  frequencies  s « s*t 
such  that 


* ii  Z,  (s)/b  • 0 {3.20} 

Shausi  and  kelly  [14]  have  generalized  the  root  locus 
techniques  of  control  theory  to  handle  characteristic 
equations  of  the  fere 
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for  distributed  parameter  syste*$  where  F5(s)  and/or  F,{s) 
can  be  written  in  an  infinite  product  expansion.  If  one  of 
the  factors  F s ( s } or  F2{s)  cannot  be  written  as  an  infinite 
product,  then  it  eust  be  a rational  function.  Equation 
(3.20)  is  of  the  for*  (3.21)  and  since  an(s)  is  expandable 
in  an  infinite  product  expansion,  we  only  require 
that  2L(s)  be  a rational  function  (i.e.,  Z,{s)  as  a “lurped 
circuit*)  or  have  an  infinite  product  expansion  'Z, (s)  in 
a distributed  parameter  system).  The  root  locus  rules 
given  by  SHausi  and  Kelly  as  modified  to  apoly  tc  the 
present  problem  of  a loaded  loop  are  given  below  [IS]. 

$TAnT  asp 

The  loci  start  ffe  * C)  at  the  zeros  of  ar(s}  and 
terminate  (K  * *}  or  the  te ros  of  2.  fs)  as  K varies  frer 
zero  to  infinity. 

Nb'XSES  3F  tDC I 

Since  **\V;  is  a trance* dental  function,  there  a ns  an 
infinite  hueber  of  loci. 

$v  tat 

The  loci  are  symmetrical  about  the  real  axis. 


The  loci  include  those  sections  o*  the  real  ***$  tHat 
lie  to  the  le*t  of  an  odd  number  oe  zeros  w*  »„  's'  arc 
l, is)  *cr  k positive  and  to  the  *t*t 


of  s r ever 


the  2eros  of  the  two  functions  for  K negative. 

iihen  a portion  of  the  real  axis  between  two  successive 
zeros  of  an(s)  or  two  successive  zeros  of  2 ^ { s ) is  part  of 
the  root  locus,  there  will  be  a particular  value  of  & for 
which  there  will  exist  a second-order  root  on  the  real 
axis.  This  root  is  known  as  the  breakaway  point. 

BREAKAWAY  P£1NT 

The  points  at  which  the  loci  break  away  fror  the  real 
axis  are  found  as  solutions  of  the  equation 


a„{s}Z,*(s)  * a‘(s)Zt(s) 


(3.22} 


ANGLES  Of  ARRIVAL  A«i&  3E9A3TU3S 

If  > ^ is  a sieple  zero  of  an(s),  the  angle  of  depar- 
ture of  tne  locus  fror  s - is  civen  by 
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If  Ss  is  i simple  sere  of  I, (s) , the  angle  of  arrival  of 


at  s*  is  giver-  by 


the  locus 


For  K > 0r  the  intersections  of  the  loci  with  the 
imaginary  axis  occur  at  those  values  for  which 

\(h)' 

arg  r = -f/2  + 2nr  n » 0,  1,  2,  . . .(3.2 

> (K). 


and  for  ii  < 0,  the  corresponding  values  of  are  found 
frcr 
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1 2 nr  n = 0,  1*2,  . . .{3. 
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S becotses  infinite*  tne  loci  ray  become  very 
complicated.  In  order  to  study  tne  loci,  we  equate  i*j 
reel  arc  imaginary  parts  of  { 3.2GJ  tc  sero: 

-**’  :»jV(*j]  - - *«  !*,{*•]  * 2 ; 3 . 
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£» inflation  of  & between;,  these  t*c  eolations  gef *•***. 
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When  Equation  (3.29)  is  used,  care  should  be  exercised 
regarding  the  sign  of  K that  has  been  eliminated.  As  written 
the  equation  define;  the  loci  for  both  positive  and  neqative 
values  of 

3.4  Step  Response  of  a Resistively  Loaded  Loop 

Figure  25  shows  the  step  response  loop  current  com- 
puted by  taking  the  Laplace  Inverse  of  (2.33)  at  } » 120° 
for  a resistively  loaded  lone  with  Zs  * * 600,  1800 

and  5400  ohms.  Figure  26  snows  the  pulse  response  at 

the  same  location  with  2 , * R,  * 600  and  1800  eh«s  for  a 
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pulse  width  equal  to  0.7 5 ctZ-'b.  Soncausal  oscillations 
can  he  seen  in  the  tine  Interval  near  t * 0 before  the 
first  signal  arrives  at  the  observation  point.  The  fre- 
quency of  the  noncausal  oscillations  corresponds  to  the 
Type  1 p©l«  of  the  first  Fourier  ssedc  nut  included  in  the 
current  representation,  as  expected.  The  various  discon- 
tinuities in  the  response  co»e  fro*  the  first  current 
Puls*  which  arrives  at  the  observation  point,  the  current 
which  travels  around  the  lonaer  path  fro*  the  source  to  the 
observation  eoint,  the  second  trip  around  the  loon,  and  so 
on.  hot*  that  for  this  loo®,  S * 13.0,  the  value  of 
ft.  * 5400  corresponds  to  a “critically  damped*  loop. 


Flo  it  re  £?> , Falser  response  atrrenf.  at  a !2Q*  for  various 

values  of  resistive  loading  far  « loop  antenna 


CHAPTER  IV 

SYNTHESIS  OF  THE  RESPONSE  OF  A LOOP  ANTENNA 


In  the  previous  chapter,  various  techniques  of  analysis 
have  been  developed  md  analytical  formulations  governing 
the  electromagnetic  behavior  of  loop  antennas  have  been 
derived.  These  techniques  and  results  are,  of  course, 

essential  in  determining  the  response  of  * loop  antenna  for 
a given  excitation  and  specifies  impedance  loading.  In 

the  following,  we  demonstrate  that  seme  ef  the  results  can 
also  be  used  to  change  the  resonant  frequencies  such  that 
a proper  excitation  {or  input),  a desired  response 
(i.e.,  current  or  radiated  field)  might  be  obtained.  Thg$ 
we  obtain,  using  the  singularity  expansion  reoresedtati on , 
a means  of  extending  to  an  electromagnetic  problem  a capa- 
bility that  is  well-developed  tr  network  theory--the  ability 
to  synthesise  the  desired  response  when  the  incut  or 
excitation  waveform  is  given. 

!rs  this  chapter  the  analytic  properties  of  the 
admittance  transfer  fuwtvlcrt  discussed  in  the  previous 
chapter  are  employed  to  formulate  the  synthesis  problem. 

5 GW®  £»f  r W»  ,•  .1*;  : ; 4 nwe  an  <4  1 - «.■»  • »**» >« 
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wentary  ti we- domain  synthesis  will  be  illustrated  by  u$inq 
the*e  findawen'al  synthesis  techniques  on  some  simple 


exaaple  problens. 

4.1  Formulation  of  the  Synthesis  Preble* 

Suppose  that  one  is  given  the  voltage  transfer*  at  the 
gap  ¥®(s),  and  one  wishes  to  find  so*e  response  quantity 
such  as  tr.c  current  at  some  point  on  the  looo  or  the  field 
at  so*e  point  in  space.  Then  for  each  Fourier  component 
of  the  current,  there  will  fee  a transfer  function  T {$} 
relating  that  Fourier  component  to  its  contribution  to  the 
response  at  the  observation  point.  The  total  observed 
response  ¥{s)  is  thew  just  the  sun 
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For  example,  if  the  current  at  a Point  ■$  is  desired,  T^Cs) 
* e^rt^t  whereas  for  a field  Point,  the  transfer  functions 
fn(si  can  fee  obtained  fro*"  Appendix  C.  Using  (3.4)  for 
the  current,  we  obtain 
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yo  *fe«) 


n*«*  in(s)  - JZL(s)/ebn 


(*.2) 


if  we  **§•*«»  that  the  input  ¥®{s)  Is  given  and  hence  fixed 
and  that  ‘tjt)  snlv  c...  i**  ***?z*^  «« 

wish  to  observe,  the  desired  response  ray  be  obtained  only 


fey  Manipulating  the  admittance  transfer  function  fey 


iepedance  loading.  Me  say  do  this  as  in  classical  network 
theory  by  prescribing  both  the  poles  and  the  residues  of 
the  edaittance  transfer  function.  The  relationship  between 
the  loading  function  and  the  poles  and  residues  has  already 
been  described  in  the  previous  chapter.  Hanety*  fro*  (3.6), 
the  condition  for  a pole  at  s 3 is 

h(*i l)  * <3-31 

and  fro»  (3.10)  the  condition  on  the  residue  requires  that 
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where  is  the  desired  residue  of  the  shifted  pole.  Thus 
the  synthesis  problem  is  that  of  requiring  interpolation 
conditions  on  Z^( s)  and  its  derivative  at  the  desired  pole 
positions.  If  the  residue  is  not  of  particular  interest, 
then  the  condition  on  the  derivative  can  be  relaxed. 

One  notes  that  this  is  not  the  usual  condition 
required  of  the  admittance  transfer  fun  r Men  in  the  cir- 
cuit theory  context . There,  the  poles  and  residues  (or 
seres)  of  the  desired  function  are  those  of  the  response 
function.  In  the  probler  posed  by  (*.3)  and  (*.*),  on 
the  Other  hand,  the  poles  and  seres  of  the  function  to  be 
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synthesized.  Z,{s),  will  not  be  those  of  the  response 
function. 

4.2  Construction  of  the  Impedance  loading  Function 
*fe  Initially  begin  with  the  simplifying  assump- 
tion that  the  impedance  loai'ng  funct ; ..  is  to  fee  synthe- 
sized using  only  linear  lusped  circuit  e lessen ts  and  devices. 
Since  we  are  uniformly  loading  the  loop,  we  way,  for  example 
use  many  electrically  small  lumped  circuits  in  series  with 
and  uniformly  spaced  around  the  loop  so  as  to  approximate 
continuous  uniform  loading.  S'ifch  the  lumped  circuit 
requirement,  Z^{s)  ®ust  be  either  a polynomial  in  $ or 
a ratio*'* ! function;  that  is,  a ratio  of  two  polynomials 
in  s.  Since  a polynomial  is  simpler,  we  consider  It  first. 

Suppose  we  have  % poles  we  wi$*  to-  synthesize  so  that 
the  condition  on  Z^($}  is  of  the  for* 

Mv)  ' " '•  J‘  • • ' . s (‘.51 

where  we  assume  that  in  the  sequence  of  soles  s,.  s2,  . . 

* If  any  pole  i$  complex.  Us  complex  conjugate  counter- 
part is  included  in  the  sequence.  The  polynomial  of 
lowest  order  satisfying  the  interpolation  condition  {4,5} 
is  constructed  using  the  Lagrange  polynomial*,  [is] 

H 
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where 


and  where 


V5>  n (> 

1 * 1 

i 1#n 


(4.7) 


(4.3) 


Thus  21  j C s .*  wiTI  he  a polynomial  of  degree  N.  - 1 , 

*f»  in  addition  to  the  pole  locations,  one  specifies 
the  residues*  we  have,  according  to  (a. 4),  additional 
constraints  on  the  derivatives  of  the  loading  function: 
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The  problem  of  interpolating  both  a function  and  its 
first  derivative  at  a set  of  points  is  solved  by  the 
Herat te  or  osculating  p olynerti  a Is  [15] 


^,s> ' I,  u-(s)  li-  * i •- 


v’ is:  z' 


In 


(4,10) 


•-here  the  functions  ( s } and  V5(*|  are  oolyno**’ als 
having  properties  similar  to  these  of  the  La  orange  inter- 


polation functions  l„(s)  of  (4.7)  and  defined  bv 
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Un{s}  « 

vn(s)  - 

tn  this  cast  Z^(s)  will  fee  of  degree  2H  - 1 , 

If  the  degree  of  the  polynomial  representing  Z,  { s) 
is  greater  than  one,  it  is  not  possible  to  synthesise 
2|Js)  using  only  passive  circuit  elements*  Since  in  »any 
applications  it  may  not  be  econcssical  !y  or  technically 
feasible  to  use  active  devices,  »e  alanine  sesre  further 
conditions  on  2^{s)  which  restrict  it  to  be  a 'posit? ve- 
rsa 1 " function  of  s . 

The  dri  vlng-point  admittance  or  impedance  functions 
of  passive  networks  {that  is*  networks  consisting  only  of 
1 taped  resistors,  capacitors*  end  inductors)  are  positive- 
real  functions*  That  Is*  our  impedance  loading  function 
m%t  be  a positive-**®*?  function  to  fee  physically  realis- 
able as  a dri vine-point  impedance  t *u*h*r  of  analytical 
properties  of  a positive-real  function  can  t*  derived 
fro*  its  definition.  The  west  basic  and  signifi- 
cant ones  are  suwaarised  in  Table  9 for  convenience, 
hole  that  {4}  in  Ttfelj  9 restricts  the  degree  of  the 
polynomial  that  can  be  used  to  represent  Z^UI  to  no 
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At  real  frequencies,  U * i**:)  the  real  part  of  tf(s)  is  an  «v*n  function  of 


greater  than  first  degree.  However,  2.  (s)  «y  be  a 

** 

rational  function  if  s. 


2l{s) 


- tP  * a -P-1 
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* a. 


fe^  ♦ b 
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*>  *-  » 


♦ b. 


{4.12} 


niter#  p * 4 * 2 « H if  Zj(s)  is  required  only  to  satisfy 
(4.3}  and  p ♦ 4 ♦ t * 2K  If  ft  also  satisfies  .M.4),  In 
either  cast,  if  JL(s)  is  to  be  positive -ret! , Ip  * qf  5,  !• 
The  coefficients  a#*  «t  * . , . , a.  ttid  be»  b$  f . . . . fer 
way  be  found  fey  substituting  (4*12?  into  (4.3 ) and  (4.4) 
and  solving  the  rt«  tilting  systee  of  linear  equations  for 
the  c-tents , 

4 ogles st ry  condition  on  the  values  2^.  that  cm  fee 
Interpolated  fey  posit? y«*raa1  functions  Has  been  devised 
fey  you!*  and  Salto  [Wife*  sail  m energy  consideration*. 

The  condition  is  t"  * the  *feeyaH«oa--r Ufc*  i U Ham* Ulan 
iitHit  Ufe»  isteMafc  emotes  yovpte*  con  legate). 


ft# 


(*  J3) 


«u$i  fee  noniiesativ*  definite,  that  is, 


* «a  * 


(4.1«! 
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for  oil  vectors  *.  Tliey  further  show  that  this  condition 
Is  Sufficient  If  the  sn  are  distinct  and  in  the  right  helf 
plane.  i.e.,Re  sn>  0,  n » 1,  2,  . . . , .1.  This  would 
Imply.  however,  that  the  loop  transfer  admittance  function 
contains  active  sources , whi ch  is  impossible.  Unfortunately , 
necessary  and  sufficient  conditions  for  the  existence  of  a 
positive  real  interpolating  function,  both  with  and  without 
the  derivative  condition,  do  not  appear  to  be  available  at 
this  time.  The  development  of  a criterion  for  the  loading 
function  to  exist  and  to  be  positive  real  is  a challenging 
problem  for  further  research. 

4.3  Time  Domain  Synthesis  Applied  to  the  Design  o*  a 

Pulse  Simulator 

One  application  of  loaded  loop  antennas  designed  to 
radiate  a specified  waveform  is  in  toe  simulation  of  the 
electromagnetic  pulse  (EMP)  generated  by  a high-altitude 
nuclear  detonation.  The  pulse  shape  required  can  be 
approximated  as  the  difference  of  two  damped  exponential 
functions,  one  having  a very  short  tie's  constant  which 
determines  the  rise  time  of  the  pulse  and  another  having 
a long  ti»e  constant  which  determines  the  rate  cr  decay 
of  the  pulse.  4 tyrical  EMP  waveform  [13]  can  be  expressed 
as 


ZA 


Best  Available  Copy 


where 

a « -2.0  x TO* 

8 * -2.6  x TO* 

The  value  of  EQ  Is  a constant,  and  for  the  purpose  of 
calcelatlons  has  been  set  to  unity. 

First*  we  wish  to  specify  the  generator  output  V*. 
For  the  pelse  generators  In  coaeon  use*  It  has  been  found 
that  generator  output  can  be  accurately  represented  by  a 
step  function  with  a finite  rise  tlae.  The  generator 
rise  tlae  and  that  of  the  wavefora  to  be  synthesized  are 
chosen  to  have  the  saae  rise  tlae  so  that  In  the  Laplace 
dona  In 


»*($)  •*{)  - e8t)  * — =£ — (*.16) 

0 ’ s(s-s) 

where  3 * -2.6  x 10*. 

In  the  far  field  there  exist  only  two  components  of 
electric  field*  E3  and  £e  (cf.  Figure  27).  The  area  of 
t»riTRa«-r  interest  for  obtaining  the  desired  transient  wave- 
* c riT  : --a  r the  axis  z*  the  loop  within  a cone  angle  of 

sbcvt  *rcr  the  axis.  The  ranges  of  angles  5 considered 
*re  the-e^ore 
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0 < 9 < 30° 
and 

ISO*  < e < ISO®  (4.171 

On  the  loop  axis  only  the  n * 1 mode  **«*ntri  bytes  to  the 
currents  and  for  the  narrow  range  of  angles  considered, 
one  may  make  the  simplifying  assumption  that  the  field 
components  are  coupled  to  only  the  n * 1 Fourier  mode  of 
current  on  the  loop.  That  is,  we  assume  that  other  modes 
of  current  on  the  loop  do  not  contribute  significantly  to 
the  far  fields  radiated  on  and  near  the  axis  of  the  loop 
antenna . 

We  wish  to  find  the  electric  and  magnetic  fields 
generated  by  a prescribed  source  function  (generator)  and 
the  resulting  induced  currents  on  ths  loop.  These  Held 
components  are  easily  obtained  with  the  help  of  the 
magnetic  vector  potential,  which  is  (cf.  Figure  27} 
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T it  complete  solution  for  the  field  components  is  obtained 
by  using  vector  di fferentis)  operators  with  (4.1£)»  An 
expansion  for  the  vector  potential  in  spherical  harmonics 
Is  derived  in  Appendix  C for  an  arbitrary  observation 
point  In  space.  This  solution  is  in  no  wsy  restricted, 

9$ 


and  near  as  well  as  far  fields  can  be  obtained.  However, 
for  the  case  being  considered,  considerable  simplification 
Is  possible  with  no  loss  in  generality  if  one  considers 
observation  points  to  be  in  the  far  field. 

The  electric  field  components  in  the  far  field  are 
given  In  [5]  for  several  nodes,  and  specifically  for 
node  n * 1, 
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where 
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for  the  range  of  e near  the  loop  axis  at  9 * 0,  (*.19)  and 
(*.20)  reduce  to 
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(4.23) 
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For  convenience  we  nil!  set  $ * 0,  factor  out  the  time 
delay  esr^c,  and  drop  all  other  unnecessary  constants, 
exhibiting  only  the  resaining  depenc^ce  of  E.  on  s: 
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Js»o(s? 

ai(s)  - jZ^s)/xbn 


(4.2*) 


Equation  (4.24)  shows  that  the  radiated  field  is  propor- 
tional to  the  ti*e  derivative  of  the  current.  For  »od« 
n • 1,  the  partial  fraction  expansion  of  the  admittance 
transfer  function  is 
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Substituting  (4.25)  into  (4.24),  we  have 
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(4.2S) 
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where  the  $j  correspond  only  to  rode  n * 1.  Thus,  it  is 
seen  that  the  effect  of  space  is  to  differentiate  the 
current  since  the  far  field  transform  is  just  proportional 
to  the  current  transform  multiplied  by  s. 

The  functional  relationship  between  the  above  quan- 
tities is  described  below: 


Adeittance  Space 

Output  * Transfer  x Transfer  x Generator 
Function  funrtion  Function 


Using  the  results  of  (4.26)  and  (4.16), 
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(4.29) 


When  the  impedance  loading  function  is  restricted  to 
be  of  the  one-  or  tw©~ele*sertt  kind*  considerable  slsrplifi 
cation  results  in  the  synthesis,  in  the  following,  atten 
tion  Is  focused  on  uniformly  distributed  resistive  and  80 
networks.  This  choice  is  cade  for  siapllcity  and  because 
resistive  and  RC  networks  are  frequently  encountered  in 
high-frequency  circuits.  Since  the  co«M>le*?ties  involved 
in  general  RLC  synthesis  are  *»uch  greater,  we  llftit 
ourselves  to  a few  basic,  si aple,  and  usefu1  techniques. 
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The  radiated  fields  of  the  unloaded  loop  antenna  do 
not  appaar  similar  to  those  of  an  EHP  waveform  due  to  the 
narked  oscillations  In  tine.  One  way  to  modify  the  radi- 
ated fields  is  to  add  resistive  loading  along  the  structure 
so  as  to  reduce  the  effects  which  cause  the  oscillations. 

If  the  structure  is  resistively  loaded  so  that  ZL($)  * 
\»  then  as  the  loading  is  increased  the  poles  in  the 
first  layer  (see  Figure  12)  wove  generally  in  the  «ob/c 
direction,  indicating  that  their  contributions  in  live 
attenuate  wore  rapidly.  The  behavior  of  the  unloaded 
Type  I pole  for  mode  n * 1,  located  where  vH/c  m 1 close 
to  the  dbfc  axis,  deserves  special  attention.  As  the 
loading  is  increased,  this  pole  moves  on  a curved  are 
down  so  the  -eb/c  axis,  at  which  point  a double  pole  is 
formed  with  its  conjugate  pole.  As  the  loading  is 
further  increased,  this  double  pole  split:-,  one  pole 
Moving  to  and  the  other  toward  zero  along  the  ob/c 
axis.  This  behavior  is  completely  analogous  to  that 
observed  as  the  resistance  is  increased  in  a series 
resonant  flic  circuit.  At  the  pof«t  where  the  double  pole 
first  is  forced,  w refer  to  v‘:«  loaded  crt  tonne  as  being 
critically  damped. 
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In  the  following,  we  examine  several  possible  approaches 
toward  the  synthesis  of  a double  exponential  waveform  (4.15). 
The  problem  might  be  considered  as  representative  of  the 
general  synthesis  problem.  In  particular,  we  encounter  cer- 
tain limitations  and  considerations  which  should  be  common 
to  any  synthesis  problem  involving  the  loop  antenna. 

The  approach  taken  here  is  to  force  the  Type  I pole  to 
be  the  synthesized  pole.  Since  it  will  have  a long  damping 
constant  and  since  the  loading  generally  forces  the  other 
poles  to  have  shorter  damping  constants  this  pole  should 
dominate  the  late  time  response.  With  the  observation  point 
along  the  loop  axis,  we  consider  only  the  n * 1 mode  and 
observe  it  in  the  far  field.  Finally  we  restrict  our  con- 
sideration to  simple  loading  functions  involving  only  resistors 
and  capacitors. 

We  begin  oy  attempting  to  specify  both  the  pole  and  its 
residue.  In  ( 4 . 1 s ) the  values  of  the  residues  are  eoual,  and 
the  requirement  exists  to  specify  only  one  remaining  pole  in 
the  su»  (4.29)  which  we  call  s’,  let  sj  be  equal  to  the 
coefficient  ot  in  (4.15),  i.e.,  the  pulse  decay  constant, 

t £ 

* -2.0  x 10  . The  remaining  task  then  is  to  equate  resi- 
dues. From  (4.29)  this  requirement  is  met  if  we  let 

3R1  g 

— • “a‘— (4.30) 

3-Sj  a^(s)  - j2L(3)/’?bn 
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where  the  impedance  loading  function  Z^(g)is  to  he  determined. 

Note  that  for  the  required  zero  of  the  transfer  imped- 
ance function,  it  is  also  true  that 


a^sj)-  jZL(s*)/^bn  « 0 . (4.31) 

Solving  equation  (4.30)  for  Z^(s)  we  obtain 
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where  is  the  value  of  the  residue  of  the  admittance  trans- 
fer function.  If  w«  specify  a series  RC  network  then 


(4.32) 


cl 


i%) 


* 


R 


±. 

Cs 


From  the  residue  condition. 
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Substituting  (4.K)  into  (4.32)  and  equating  to  (4.33) 
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Similarly,  substituting  (4.33)  Into  (4.3!)  yields 


VM> 


(4.36) 


Thus  we  hive  two  equations  with  two  unknowns,  which  can  be 
solved  for  R and  C. 

Solving  the  system  of  equations  (4.3S)  and  (4,36}  gives 
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1 


j*bn  [(a  - sj)a'(sj)  ♦ a(sj)  - a(8)] 


(4.37} 


and 


R * (s^ ) 


(4.38) 


The  solution  expressed  by  (4.37)  and  (4.3i#)  is  theoretically 
correct;  in  practice  it  is  nr  realisable  with  passive  ele- 
ments since  the  denominator  of  C is  negative,  whereas  the 
numerator  is  positive  Thus  specifying  Soth  the  pole  and  its 
residue  yields  an  unphysical  solution. 

One  sees  that  the  requirement  that  the  residues  be  equal 
arises  from  (4.?S)  because  the  response  at  t * 0 should  ae 
aero,  Since  the  short  time  constant  exponential  is  provided 
by  the  source  and  the  longer  time  constant  imes  fror  the 


antenna  "ringing  down,"  and  the  wavefom  shape  during  the 
transition  between  the  rise  and  decay  is  not  critical,  one 
should  be  able  to  obtain  roughly  the  desired  response  with- 
out specifyinq  the  residue  of  the  pole.  Accordingly,  three 
additional  cases  were  selected  where  the  pole  location  was 
specified  to  yield  a decay  constant  equal  to  the  value  of  a 
soecified  in  (4.H).  The  first  case  used  a ourely  resistive 
load  of  4944  H;  in  the  second  case,  3000  r.  and  66  uf  capac- 
itance were  used;  and,  in  the  third  case,  a 1000  a resis- 
tance and  399  uf  capacitance  were  used.  These  cowbinatiors 
were  chosen  by  requiring  2^Uj)  * 4944  ♦ 10  which  puts  the 
pole  at  the  desired  position.  The  ti«e  domain  response  for 
each  case  is  plotted  in  figure  28.  In  the  figure,  the  shape 
in  each  case  is  almost  ident'cal,  as  expected.  However, 
the  decay  tin*  is  ouch  shorter  than  the  desired  value.  Tp 
see  that  this  effect  is  Independent  of  both  the  generator 
excitation  pulse  shapr  and  the  loading,  the  step  function 
response  was  compute4  and  is  shown  in  Moure  2?.  Tht  simi- 
larity of  the  response  in  each  case  leads  one  to  the  con- 
clusion that  it  is  a sere  of  the  loop  transfer  ad*M  ttance 
it  i • 0 which  causes  the  difficulty.  Recall  that  the 
unloaded  loot?  transfer  impedance  function  has  a pole  *t 
$ * 9 which  translates  to  a cere  i«  the  loco  admittance  for 
both  the  loaded  and  the  unloaded  case.  This  transmission 

zero  tends  to  cancel  the  synthesized  pole  at  $ » «: 
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’'aoeda^te  Ic^^ps. 


which  Is  c lose  tc  the  origin  resulting  in  a very  decreased 
late  tint  response. 

This  rero  arises,  of  course,  from  the  combined  ef- 
fect ef  the  remaining  poles  of  the  loop  which  have  been 
unconstrained.  Hence,  thefr  effect  in  the  late  t i see  re- 
sponse is  not  negligible  as  was  originally  assumed. 

Since  a^fs)  has  a pole  wt  s * g then  for  low  fre- 
quencies the  loop  transfer  impedance  is  capacitive.  Since 
the  unloaded  loop  is  passive,  this  equivalent  circuit  ele- 
ment must  be  positive.  In  orinr  for  loading  to  be  chose* 
so  as  to  came*,  the  pole  of  a. is)  (i.e.,  the  tero  of  the 
admittance  transfer  function),  l.ii)  would  have  to  cancel 

1m 

the  low  frequency  behavior  of  a U}*  which  would  reoulre  a 
non-physical  rtgative  capacitance.  Hence,  it  is  rot 
possible  to  cancel  the  sees  in  the  admittance  transfer 
function  b?  us f fig  passive  leading. 

Td  test  the  validity  of  this  «srlan*t«oa  for  the 
poor  site  time  behavior,  a numerical  experiment  was  con- 
ducted id  determine  If  H was  possible  to  eliminate 
the  **re  in  the  admittance  transfer  function,  accord- 
ingly the  generator  output  waveform  was  modified  to  the 
time  integral  of  the  original  excitation,  which  intro- 
duced another  factor  of  Vs  in  the  transform  domain  so 
is  te  cancel  the  tare  in  the  transfer  admittance  at  s « V 
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This,  of  course.  Introduces  « reap  in  the  tlae  doaeln 
output  of  the  generator  but  fulfills  our  reoulreaent 
for  an  additional  pole  In  the  denoalnator  of  (4.27). 

The  resulting  tlae  domain  response  is  plotted  in  Figure  30. 
The  similarity  of  the  late  tlae  response  to  the  desired 
double  exponential  conflraed  the  conclusion  that  the  zero 
in  the  admittance  transfer  function  caused  the  previous 
difficulty  in  achieving  good  late  time  response.  The 
integration,  however,  further  degrades  the  early  tlae 
response. 

Since  good  early  tlae  response  was  obtained  with  a 
generator  frequency  dependence  given  by 

V*(s)  , — 2*_ 

0 s(s-3)  (a. 39) 

while  good  late  tine  doaair.  behavior  was  obtained  with 
the  '-•suercy  dependence 


sMs-3) 


(4.40) 


one  night  speculate  that  perhaps  a good  overall  approxi- 
mation tc  the  desired  ^esaonse  night  be  obtained  by  the 

p V - ^ ^ n r, 
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• A Vnlf^ui 


'if  'ie’s  far 

aoai f 1««5  ra.r»o  far  various 

values  of  laattfaoci  ioi«H«s. 


(♦.*1) 


ve(s)  * ziilDl 
0 s2(s-8) 

which  Has  the  desired  high  frequency  (early  tine)  behavior 
of  (4.39)  and  the  desired  low  frequency  (late  time)  be- 
havior of  (4.40).  That  this  Is  Indeed  the  case  Is  seen 
by  noting  that  the  response  of  (4.41)  Is  obtained  by  super- 
position of  the  responses  due  to  (4.39)  (figure  28)  and 
(4.40)  (Figure  30).  However,  this  approach  (i.e.,  modify- 
Ing  the  generator  waveform  to  produce  the  desired  response) 
runs  counter  to  the  objective  of  synthesizing  the  desired 
response  by  loading  the  loot.  It  appears,  then,  that 
although  we  way  be  able  to  synthesize  the  pole  pattern  of 
the  loop  for  a finite  number  of  the  poles,  we  eay  require 
a wore  elaborate  treatment  to  guarantee  that  the  petition- 
ing of  a finite  nuwber  of  the  poles  by  Impedance  loading 
does  indeed  lead  to  the  desired  time  domain  waveform. 

The  complexities  introduced  by  the  infinite  number 
of  poles  and  the  apparent  late  time  differentiation  of 
the  decaying  waveform  due  to  a zero  in  the  admittance  trans- 
fer function  pose  a unique  and  difficult  set  of  constraints. 
This  is  a problem  outside  the  scope  of  this  study  but 
offers  interesting  possibilities  for  future  research. 
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CHAPTER  V 
CONCLUSIONS 

The  objective  of  this  research  has  been  to  develop  some 
fundamental  techniques  for  the  analysis  and  synthesis  of  the 
response  of  a loaded  loop  antenna.  In  the  past,  the  tine  do- 
main raspon'u  for  such  a problem  would  be  determined  either 
by  time  harmonic  analysis  coupled  with  Fourier  inversion, 
or  by  direct  time  domain  solution.  Hlth  the  addition  of  im- 
pedance load tug,  considerable  effort  would  then  be  spent  re- 
calculating the  entire  response  of  the  loop  antenna  without 
making  use  of  any  of  the  information  about  the  response  of 
the  unloaded  antenna.  Use  of  the  singularity  expansion  Method 
(SEMI,  however,  permits  one  to  systematically  examine  the 

effects  of  loading  using  the  solution  for  the  unloaded  loop 
antenna . 

The  observation  that  the  solutions  of  electromagnetic 
problems  are  analytic  functions  of  the  complex  frequency  s 
except  at  singularities  forms  the  basis  of  the  SEW  and  per- 
mits one  to  use  the  many  powerful  theorems  of  complex  varia- 
bles to  sore  efficient’'  represent  the  solution.  The  result- 
ing time  domain  response  representation  is  a superposition 
of  damped  exponentials  whose  complex  freouencies  correspond 
to  s-pUne  poles  of  the  admittance  transfer  function.  These 


ZZ'ZZ  ««tr  deieranneti  »r*na i fcne  impedance  loading  and  the  es~ 
leaded  admittance  transfer  function.  Thus,  the  sdvaetage 
of  the  singularity  expansion  technique  is  that  one  can  sepa- 
rate and  characterize  basic  attributes  of  the  structure  only 
once*  and  the  tine  doaiain  response  for  various  loadings  and 
excitations  can  then  be  easily  determined  freat  the  structure’s 
characteristic  behavior. 

In  Chapter  I 15*  it  is  shown  that  s-plane  contour  plots 
of  the  magnitude  and  phase  of  the  unloaded  impedance  transfer 
function  of  the  loop  permit  one  to  readily  determine  the  tra- 
jectories of  the  poles  as  loading  1$  added  to  the  structure. 
Furthermore,  the  observation  that  the  loading  can  be  inter- 
preted os  adding  a feedback  path  to  the  admittance  transfer 
function  permits  one  to  use  the  root  locus  techniques  of 
control  systems  to  further  aid  in  the  determination  of  the 
pole  movements  with  increased  loading.  Since  we  are  dealing 
with  an  antenna  that  is  a distributed  parameter  system,  the 
conventional  root  locus  technique  was  generalized  so  as  to 
be  applicable  to  « system  with  a countably  infinite  number 
of  singularities.  The  generalized  root-locus  technique  pro- 
vides a valuable  tool  for  studying  the  effect  of  varying  the 
impedance  loading  over  a wide  range. 

The  synthesis  of  time  domain  waveforms  by  impedance 
loading  has  been  considered  in  Chapter  IV.  It  was  found  that 
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the  required  condition  on  the  impedance  loading  to  locate 
a pole  at  any  point  in  the  complex  s-plane  is  that  it  Inter- 
polate the  impedance  transfer  function  at  the  desired  pole 
frequency.  If  the  residue  at  the  ode  is  a1  so  to  be  speci- 
fied, the  derivative  of  the  impedance  function  also  satisfies 
an  interpolator  constraint.  These  conditions  may  be  satis- 
fied by  Lagrange  or  Hermite  interpolating  polynomials,  res- 
pectively. However,  if  one  is  restricted  to  passive  loading, 
the  loading  function  must  be  a rational  function  of  s.  A 
necessary  condition  on  the  interpolation  constraints  is  given 
for  the  realization  of  passive  loading.  A sufficiency  con- 
dition for  realization  with  passive  elements  is  apparently 
lacking  at  this  time  However.  Some  sime1*  cU&mpts  in  Chap- 
ter IV  to  synthesize  a radiated  waveform  consisting  of  the 
sue  of  two  exponential  functions  were  only  partially  succ^  s- 
ful.  The  difficulties  seemed  to  arise  froa  attempting  to 
control  an  infinite  number  of  poles  by  loading  and  from  the 
presence  of  zeros  in  tne  output  response  due  to  the  admittance 
transfer  function  and  the  free  space  transfer  function. 

In  summary,  this  research  was  directed  towards  simpli- 
fying the  understanding  of  impedance  loaded  lot,*  antennas 
using  the  singularity  expansion  solution  technique.  On  the 
basis  of  the  results  of  this  study,  several  recommendations 
concerning  future  research  are  suggested.  Furthe/  study 

:?5 


1$  needed  to  determine  whet  constraints  exist  on  the 
realizability  of  passive  loads.  Also  needed  is  an  approx- 
1«ate  theory  for  treating  the  infinite  nuaber  of  poles 
of  Type  III.  Possibly  a transmission  line  model  would 
enable  one  to  factor  these  poles  out  of  the  admittance 
transfer  function  enabling  one  to  work  w’'th  a few  oolet  in 
partial  fraction  for*,  the  rest  being  absorbed  into  a 
transcendental  function  representing  the  transmission  line 
approximation.  Hopefully,  such  an  approach  might  lead  te 
a better  understanding  of  the  constraints  on  the  realizabil- 
ity of  time  domain  waveforms  Imposed  by  the  structure. 
Another  approach  to  .synthesis  might  involve  optiaization 
techniques  to  choose  the  loading  so  as  to  minimize  the  error 
between  the  desired  response  and  that  actually  obtained  from 
the  antenna.  Here  the  shifted  poles  would  iso t be  specified 
but  would  enter  the  calculations  only  as  a means  to  compute 
the  time  domain  response.  Finally,  an  interesting  area 
for  future  research  is  in  developing  efficient  ways  to  han- 
dle non-uni  form  and  point  loading.  In  the  case  of  the  loop, 
such  loading  unfortunately  couples  all  the  modes  together. 
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APPENDIX  A 

DERIVATION  OF  THE  INFINITE  PRODUCT 
REPRESENTATION 

tteierstrass'  Theorem  for  infinite  products  [19]  requires 
derivatives  with  respect  to  s.  In  the  following,  it  is  nota- 
tionally  convenient  to  use  both  the  wavenumber  k » *j$/c  and 
the  Lapi  re  transform  variable  s simultaneously.  Thus,  for 
example, 

*»*<«> . *a.-i  li. 

ds  dk  ds  c dk 

Except  for  n « 0,  *M(s)  has  a pole  at  s « 0 which  we  wish  to 
eliminate.  Hence  we  consider  the  intermediate  function 

ffi(s)  * saft(s),  n * *1,  +2,  . . . C A~ 2 ) 

which  has  only  zeros  in  the  finite  complex  plane  and  hence  is 
an  entire  function. 

lie  now  consider  the  logarithmic  derivative 

f;U)  # aR{$)  ♦ s*;(s) 

f„(s)  sa^is)  (A-3) 

which  is  meromorphic.  If  this  function  is  bounded  on  a set 
of  contours  Cn  enclosing  the  poles,  then  an  infinite  product 
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The  function  F ( z ) Is  defined  as 

/2  j 

[dn{2)  - jfinU)]dz 


and  the  constant  C„  is 
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C * In  4 ♦ y-2]T  L~». 

n?o  > 


(a-b; 


(a-t: 


U-S) 


In  (A-6),  hQ  and  I0  are  Modified  Bessel  functions,  y is 

Euler's  constant,  and  z * ‘iks 

£ 

The  function  in  (A-7)  Fw(z)  can  also  b®  written  in  integral 
fora  as  [21] 
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Next  we  derive  a recursion  fonaula  for  the  derivative 


119 


j[z  sin©  - (n-1)0] 


Fa{2).  Froa  (A«9),  we  have 
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Integrating  the  second  ter*  yields  finally 
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Differentiating  ( A- 9 ) with  respect  to  z,  we  have 


F'(z)  « 1 
n * 


- j C 2 Sine  -r»6  } 

de 


which  upon  comparison  with  C A- 1 6 ) yields  the  desired 
formula. 
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Using  the  Fundamental  Theorem  of  Calculus  with  (A-7) 
alternatively  write 
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where  U-18)  has  been  used  In  (A-23). 

Uslns  the  asyaptotic  expansion  obtained  by  Uaashaiikar  [8J 
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so  that  finally,  wo  have 
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Therefore  {A- 3)  reduces  to 
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TM$  function  is  bounded  in  both  the  left  and  right  $ -plan® 
as  s~*  and  also  on  a circular  contour  which  passei  between 
the  poles  of  f „ ( ;s } which  are  the  zeros  of  *n(s). 

Th?>$,  the  product  expansion  of  ?*U)  is  given  by  [20] 
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Finally,  we  have  that 
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so  that  eoabining  (A-33)  - (A- 39 ) . we  have 
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Thus,  the  infinite  product  representation  for  an C 4 > is 
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Turning  to  the  representation  of  a0(s),  we  note  tha& 
a0(s)  has  a sero  at  s 0,  so  we  consider  the  function 
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The  asyaptotic  formulas  froa  [8],  yield 
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This  function  Is  bounded  on  a circle  of  radius  R attached  at 


theoriginand  passing  between  poles  In  the  asymptotic  layer. 
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Hence  — — - Is  bounded  on  s sequence  of  such  circles  (L.  en- 
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closing  p poles.  He  now  need  the  value  of 


fi<0) 


*;(s)  . i\ 
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Fro*  (A-55)»  we  have 
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so  that 


^0)  . Q 

fjo) 
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and  we  have  finally* 
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where  the  produce  Is  over  all  the  zeros  of  &ni%)  except  s*0. 


Since 


lin  «0(s) 

fo(0)  “ 4fcf 


0(s)  1 1b  k&K1 
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we  ^ave  finally. 
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CALCULATION  OF  THE  NATURAL  FREQUENCIES  BY  TNE  METHOD 
OF  MOMENTS  FOR  THE  n * 0 HCOE 


The  tsethed  of  sscsaenfcs  solution  for  a ode  * 0 requires 
oely  a component  of  surface  current  which  is  $ in defend 
ent.  By  sysoetry , there  exists  only  & $ c exponent  of 
eagnetic  vector  potential  which  is  also  $ independent.  Th 
scattered  electric  field  is  given  by 


j«U0£ 


(k*  ♦ ??  • ) A 


(B-1) 


P 
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£ - <&; 


i < \jj 


Because  of  the  a Independence*  the  a cowponent  of  the 


scattered  electric  field  is 

t m A 

0 j*U£  * 

where  the  Baguette  vector  potential  is  given 


(B-2) 
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y p*  -IjcR 
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and  where 


o'dd'adt'1 


o o 


(6*3) 


-it  vs'!*  , r'1* S'  4'v 


R1  • |r  - ?' | 1 


2a8(l  - sin  $ sin  $a  - cos  $ cos  cos  o') 


♦ 2b*(l  - cos  o') 


♦ sJabO  - cos  $‘)(cos  o * cos  $’)  (6-4) 


o'  * b ♦ a cos  v ' 

The  coordinates  o'*  and  $'  are  define^  In  Figure  1. 
Since  the  fields  are  o Independent,  e Has  been  set  eaual  to 
rero  In  (8-3)  and  (8-4).  The  singularity  occurring  in 
(B-3)  when  o * O'  and  c ’ • 0 (1  .e. , R • 0)  Is  difficult  to 
Handle  In  a numerical  solution.  Accordingly,  we  attract 
the  singular  part  of  the  integrand  analytically  In  the 
following. 

Considering  the  Integration  on  o'  first,  we  write 
the  distance  from  the  source  point  to  the  field  point  as 


R « (B  - C cos 


(8-S) 


whare 


8 « 2a1  - Ea*  sin  o 


O'  * 2b2  * 2«b(co$  o ♦ cos  #*) 


C * Sa2  cos  o cos  v'  * 2b3  + 2*b(cos  « * cos  */')  (6-6 ) 
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so  that  (B-3)  becooes 


ZS 

!a  f 

2*  1 


«J^($')(b  ♦ a cos  $')  a 


-jk(8-C  cos  o,)i/* 
cos  $*  S do*  <H* 

(B  - C cos  y),/a 


Consider  the  Integral  in  the  brackets.  Ut  isolate  the 
singularity  by  adding  and  subtracting  a ter*  having  the 
sa*e  singularity  as  the  Integrand  hut  which  is  Integrate 


cos  ($*) 


• jk(8  - C cos  $*P 
(9  - C cos  d1®* 


MB-C  cos  *‘}‘/S 
(8  - C cos 


* I cos  U ) 2 — -»-• ■ — — L 


:«*).  dd* 


(B-8) 


% (B  - C cos  a*)  '* 


The  first  integral  on  the  right-hand  side  is  nonsingular 


m 


and  is  hence  aeenable  to  nuatericel  integration.  The  second 
integral,  which  we  now  proceec  to  evaluate,  contains  the 
slngulari ty. 

Mith  the  substitution  cos  • (Sees2  o'/?  - ’)  and 
the  change  of  variables  $*  * * - 2C,  the  second  integral 
bee owes 


CCS  o'  do" 

(3  - C COS  o' )‘/S 

»/2 

. 2 I (1  ~ n sift*  U *■  C«/2  - VI 

{3  * 4 -a/2  (1  - s»  sin2  S)1^* 
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r / _ v <L 
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( 1 - 2/*) 
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(»(B  * C) 


(8  ♦ C)l/S 


where  w * 2C/B+C  and  K(~)  and  £(r)  are  elliptic  integrals 
of  the  first  and  second  kind,  respectively  [22],  Thus, 
with  (S-8)  and  (8-9).  (8-7)  becowes 


«Mv‘ ) (b  * a cos  *’  5 a 
o 


-JM9-C  cos 


(s  - c cos  y)‘ 


«(8  ♦ C) 


CS  ♦ £) 


The  tern  involving  K(«)  Is  still  singular*  since  as  t 
approaches  v’ . ■ tends  to  unity.  Using  (8 -6},  «ie  rewrite 
is  to  exhibit  its  uependenee  on  *;•'  explicitly; 


D » F cos  v* 

F * G COS  t ' - H Sin  it ' 


(3-11 ) 


where 


0 * 4hs  ♦ 4ab  cos  v 


t * *as  cos  v * Sab 


¥ * 2&3  * 4frl  ♦ 4«<j  ccs 


0 * 2a1  cos  v * 4«o 


■■ 


H « 2&*  sin  o- 


The  singular  integral  of  interest  is 


0 Jy'  ) {b  ♦ a cos  % * ) a 


>2(1  - 2f S 
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K(«)  d* 


, ,:-A>  s’. 


(5-13) 


V*,,(b * a cos 
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in  » fc  f «»  j 

(B  ♦ 0,/S 


- (-2  In  j $ - $•  |)  { B - 1 6 ) 


Vie  now  consider  dividing  the  cross  section  of  the  loop  Into 
N subsections  of  angular  extent 


and  define  eidpcints  and  end  points  o*  each  interval  as 


(n  - 1)  i; 


vn«  ■ (n  - 1/2 ) Cv 


2n.  “ (n  - 3/2}  iv  n a 1,  2, 


The  current  is  expanded  in  pulse  functions 


. . 1 


(8-18) 
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where 
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$ Is  3 


1 * &*-  < $ <.  $„♦ 


P,t»! 


0,  otherwise 


and  substituted  la to  the  vector  potential.  Since  by  (B-2), 


vector  potential  Is  proportional  to  the  electric  field,  at 


a natural  resonant  frequency,  the  vector  potential  due  to 


the  current  along  the  surface  1$  zero.  If  this  condition 


Is  enforced  at  the  points  vp . p ■ 1,  2,  . . . , N,  a eatrlx 


results  whose  deterwinant  Is  zero  at  the  pole  frequency. 


That  is, 


det  !l(s) ! * 0 


(*-21) 


when  s Is  a natural  resonant  frequency.  The  natrlx  Z(s) 


1$  defined  by 
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(3  • C cos  o’)1/J 
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2 In  (bp  - ^j)  d*' 
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where  Mb')  » '.  ♦ * cos  p'  in  the  above  expressions,  “ofce 


that  the  singular  tore  (5-16)  has  been  extracted  fron  the 


Integrand  and  its  integral  Is  added  outside  the  integral  in 


(8-22)  for  p “ n, 


The  Fourier  expansion  of  the  voltage  across  a uni  for* 


gap  Is 
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Using  (8-22) , we  can  write  a function  dependence  for  the 


total  current  where 


JLL  Z I - V 
juue 


(B-24) 


It  Is  pointed  out  that  the  total  current  Is  equal  to  th® 


sue  of  Individual  current.  That  Is, 


Lln  • e 


•o<*> 


B-25) 


where  C * -uV/tn^  and  the  sun  Is  the  quantity  plotted  in 
Chapter  II,  Figure  10. 
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"vsAilWBk.^  (5 
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nay  be  determined  free  the  vector  potential  for  any 
arbitrary  point. 

The  element  of  vector  potential  d A at  a point  r 
0 or  S,  Y,  T,  figure  27,  has  two  components 

dAx  » -dA  sin  o' 

dAy  * dJI  cos  o'  ( 

Those  eay  he  expressed  as 

a*  . ~~£.  i (e* ) S sin  o'  bdo* 

* 4*  R 


w -jkR 

a*.,  - !(♦')  £—  cos  e*  bdo*  ( 

* 4*  8 

He  tay  writ** 
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where  coops  ring  with  ( C-1 ) , we  have 
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Furthcrsore,  noting  that 
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5 - -Jkh<Z,(kR) 

R 0 
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where  hj  Is  the  spherical  Hanfce'  function  of  order 
zero,  second  kind,  we  ray  employ  the  addition  theorem  [23] 
to  write 


hJ(kR)  « £ (2t*l)  h[2*(kb)  J,(krQ)  P£(cos  C)6  rQ  < b 
t *®  0 


£ (21+1)  h[£^(kr0)  J£{k&)  P£{cos  U . rQ  > h 

cc-n 


The  Legendre  functions  can  further  be  expanded  as  [23] 


P.icos  U . t e.  P'fcos  0! 

£ «M)  0 (£♦«)!  1 


P"{ces  e’)  cos  «Ko-c*}  (C-8) 
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ufetrt  c#  * ! for  * • 0 and  ea  * 2 for  a > 0.  This  Bay 
also  be  written  as 


’.(cm  o ■ £ IkialiL  p!«l(cos  e) 

£-1  (t*|«|)!  1 


p|b! (cos  6')  *-■*■<♦-♦')  (C-9) 


Hence,  for  re  > is',  we  have 
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Substituting  (€-9)  Into  (C-10),  we  way  rearrange  the  order 
of  the  suwaation  to  obtain 
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Let  the  inner  suaaatfnn  be  represented  by  an  indexed  tern 


J e ( kb } 
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( C-12 } 


so  that  the  vector  potential  can  be  written  as 
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Integrating  over  all  o’  and  using  the  orthogonal i ty  of  the 
function  @xp(j«-d'},  w*  obtsin  finally. 
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